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Abstract. Parametrization of 4 x 4-matriccs G of the complex linear group GL(4,C) in 
terms of four complex 4- vector parameters (k, m, n, I) is investigated. Additional restrictions 
separating some subgroups of GL(4, C) are given explicitly. In the given parametrization, the 
problem of inverting any 4x4 matrix G is solved. Expression for determinant of any matrix G 
is found: dct G = F(k, m, n, I). Unitarity conditions G + = G _1 have been formulated in the 
form of non-linear cubic algebraic equations including complex conjugation. Several simplest 
solutions of these unitarity equations have been found: three 2-parametric subgroups G\, 
Gii G3 - each of subgroups consists of two commuting Abelian unitary groups; 4-parametric 
unitary subgroup consisting of a product of a 3-parametric group isomorphic SU(2) and 1- 
paramctric Abelian group. The Dirac basis of generators A^, being of Gell-Mann type, 
substantially differs from the basis used in the literature on SU{4) group, formulas 
relating them are found - they permit to separate SU (3) subgroup in SU (4). Special way to 
list 15 Dirac generators of GL{A, C) can be used {A k } = {cvj©/3j © {cuVPj = K®L®M)}, 
which permit to factorize SU{4) transformations according to S = e laa e lb P e lkK e llL e lmM , 
where two first factors commute with each other and are isomorphic to SU{2) group, the 
three last ones are 3-parametric groups, each of them consisting of three Abelian commuting 
unitary subgroups. Besides, the structure of fifteen Dirac matrices A& permits to separate 
twenty 3-parametric subgroups in SU{4) isomorphic to SU{2); those subgroups might be 
used as bigger elementary blocks in constructing of a general transformation SU(A). It 
is shown how one can specify the present approach for the pseudounitary group SU(2, 2) 
and SU(3,1). 
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1 Introduction 

The unitary groups play an important role in numerous research areas: quantum theory of 
massless particles, cosmology models, quantum systems with dynamical symmetry, nano-scale 
physics, numerical calculations concerning entanglement and other quantum information pa- 
rameters, high-energy particle theory - let us just specify these several points: 

• SU {2,2) and conformal symmetry, massless particles [71 [T9 l [2tj j [39 | 172]; 

• classical Yang-Mills equations and gauge fields |64j ; 

• quantum computation and control, density matrices for entangled states [2| I31|. [65] : 

• geometric phases and invariants for multi- level quantum systems |55j ; 

*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
|http: / /www.emis.de/journals / SIGMA /symmetry2007.html | 
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• high-temperature superconductivity and antiferromagnets [361 152]; 

• composite structure of quarks and leptons [67 1 f68l I5T] ; 

• 5*7(4) gauge models [731 [29]; 

• classification of hadrons and their interactions [30 1 134 1 [25]. 

Because of so many applications in physics, various parametrizations for the group elements 
of unitary group 5'C/(4) and related to it deserve special attention. Our efforts will be given to 
extending some classical technical approaches proving their effectiveness in simple cases of the 
linear and unitary groups SL(2, C) and SU(2), so that we will work with objects known by every 
physicist, such as Pauli and Dirac matrices. This paper, written for physicists, is self-contained 
in that it does not require any previous knowledge of the subject nor any advanced mathematics. 

Let us start with the known example of spinor covering for complex Lorentz group: consider 
the 8-parametric 4x4 matrices in the quasi diagonal form [181 

ko + ka 

G = 

?tt,q — ma 

The composition rules for parameters k = (ko,k) and m = (mo,m) are 

/lq /uq/cq - \~ k k ^ k Jxi^k ~\~ k JvQ ~\~ %k x k . 

m,Q = m' mo + mm, rn" = m' Q m + rn'm^ — im' x rn. 

With two additional constraints on 8 quantities k$ — k 2 = +1, ?h,q — m 2 = +1, we will arrive at 
a definite way to parameterize a double (spinor) covering for complex Lorentz group 5'0(4, C). 
At this, the problem of inverting of the G matrices with unit determinant det G is solved straight- 
forwardly: G = G(ko,k,mo,m), G" 1 = G(ko, — k, uiq, — m). Transition from covering 4-spinor 
transformations to 4- vector ones is performed through the known relationship Gj a G~ 1 = ^ c L c a 
which determine 2 =^ 1 map from ±G to L. 

There exists a direct connection between the above 4-dimensional vector parametrization 
of the spinor group G(k a ,m a ) and the Fedorov's parametrization [32] of the group of complex 
orthogonal Lorentz transformations in terms of 3-dimensional vectors Q = k/ko, M = m/mo, 
with the simple composition rules for vector parameters 

„ _ Q + Q' + iQ' xQ , _ M + M' — iM' x M 

~ l + Q'Q ' ~ 1 + M'M 

Evidently, the pair (Q, M) provides us with possibility to parameterize correctly orthogonal 
matrices only. Instead, the (k a ,m a ) represent correct parameters for the spinor covering group. 
When we are interested only in local properties of the spinor representations, no substantial 
differences between orthogonal groups and their spinor coverings exist. However, in opposite 
cases global difference between orthogonal and spinor groups may be very substantial as well as 
correct parametrization of them. 

Restrictions specifying the spinor coverings for orthogonal subgroups are well known [32] . 
In particular, restriction to real Lorentz group 0(3, 1) is achieved by imposing one condition 
(including complex conjugation) (k,m) ==^ (k,k*). The case of real orthogonal group 0(4) is 
achieved by a formal change (transition to real parameters) (ko,k) ==^ (ko,ik), (mo,m) =^ 
(mo,im), and the real orthogonal group 0(2, 2) corresponds to transition to real parameters 
according to (k ,ki,k 2 ,k 3 ) ==> (k ,ki,k 2 ,ik 3 ), (m ,mi,m 2 ,m 3 ) ==> (m , mi, m 2 , im 3 ). 

To parameterize 4-spinor and 4-vector transformations of the complex Lorentz group one 
may use curvilinear coordinates. The simplest and widely used ones are Euler's complex angles 
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(see |32] and references in |18j). In general, on the basis of the analysis given by Olevskiy [58] 
about coordinates in the real Lobachevski space, one can propose 34 different complex coordinate 
systems appropriate to parameterize the complex Lorentz group and its double covering. 

A particular, Euler angle parametrization is closely connected with cylindrical coordinates 
on the complex 3-sphere, one of 34 possible coordinates. Such complex cylindrical coordinates 
can be introduced by the following relations |18j : 



k$ = cos pcos z, &3 = i cos p sin z, k\ = i sin p cos (f), k 2 = i sinpsin^), 
mo = cos R cos Z, 7773 = i cos R sin Z, m\ = i sin 777,2 = i sin R sin <&. 

Here 6 complex variables are independent, (p, z, <$>), (R, Z, $), additional restrictions are satisfied 
identically by definition. Instead of cylindrical coordinates in (p, z, <f>) and (R, Z, <!>) one can 
introduce Euler's complex variables (a,f3, 7) and (A,B,T) through the simple linear formulas: 



a 



+ z, P = 2p, j = (j)-z, A = $> + Z, B = 2R, 



<f>-Z. 



Euler's angles (a, (3, 7) and (A,B,T) are referred to k a , m a -parameters by the formulas (see 
in [32]) 



cos (3 - 
cos a ■- 

cos 7 = 

cos B 
cos A 

cos r = 



— k^ + k\ + k^, 



-ik ki + k 2 k 3 



V^o V~ ^1 ^2 
-ik ki - k 2 k 3 



k 2 



2 2 1 2 1 2 

TTIq — 777g + 771]^ + 7Ti 2 , 

+«momi + m 2 m 3 



sin/3 



+imQmi — 772-27713 



llln 



777: 



IV 1 



■777T 



777S 



A- 2 
«, 2 , 



sin a 
sin 7 : 
sin 73 : 



v 7 ^ 



k 2 



2 

sin ^4 
sinT : 



7?)? 



772? 



-7277 



777. 



2' 



+777707772 — 777i?773 



/?)? 



777: 



-7777 



777? 



+777707772 + 777i7773 



777? 



777: 



7777 



777? 



Complex Euler's angles as parameters for complex Lorentz group <SO(4, C) have a distin- 
guished feature: 2-spinor constituents are factorized into three elementary Euler's transforms 
{a 1 stands for the known Pauli matrices): 



B(k)-- 
5(777) 



e -i^a/2 e ia^/2 e +ia^/2 g S L(2,C), 

-- e — 3 r/2^B/2 e WA/2 e SL (2,C)'. 



The main question is how to extend possible parameterizations of small orthogonal group 
50(4,(7) and its double covering to bigger orthogonal and unitary groupd3- To be concrete we 
are going to focus attention mainly on the group SU(4) and its counterparts SU(2, 2), 5(7(3, 1). 

There exist many publications on the subject, a great deal of facts are known - in the following 
we will be turning to them. A good classification of different approaches in parameterizing finite 
transformations of ST/ (4) was done in the recent paper by A. Gsponer [35] • Recalling it, we will 
try to cite publications in appropriate places though many of them should be placed in several 
different subclasses - it is natural because all approaches are closely connected to each other. 



1 In this subject, especially concerned with generalized Euler angles, we have found out much from Murnaghan's 
book [561 . 
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• Canonical form [531 El E3 E7J EH El [75]. They use explicitly the full set of the 
Lie generators^] so that the group element is expressed as the exponential of the linear 
combination 

G = exp n [i(aiAi H h a n \ N )} 

the infinite series of terms implied be exp - symbol is usually very difficult to be summed 
in closed form - though there exists many interesting examples of those: 

• Non-canonical forms [3l[l[g[9l[IIl[3g[38j[^ As a consequence 
of the Baker-Campbell-Hausdorff theorem [H [26], [37] it is possible to break-down the 
canonical form into a product 

G = exp n( i) x • • • x exp n(fe) , H h n {k) = N. 

with the hope that exp n (i) could be summed in closed form and also that these factors have 
simple properties. This possibility for the groups S'f7(4) and SU {2,2) will be discussed in 
more detail in sections below. 

• Product form |12[ [13] [T4"l [T6l [271 l38l l56l 157] . An extreme non-canonical form is to 
factorize the general exponential into a product of n simplest 1-parametric exponentials 

G = expfiaiAi] x • • • x exp[iaArAjv]. 

• Basic elements (the main approach in the present treatment) [HI [H [351 HOI EH H3 EH 
Wf\ 174] . This way is to expand the elements of the group (matrices or quaternions) into 
a sum over basis elements and to work with a linear decomposition of the matrices over 
basic ones: 

G' = x' n X m , G = x n X m , X = I, k G {0,1,... ,N}, 

G — G G, XfcXfc — x m A m x n A n — x m x n \ m X n , (!•!) 

as by definition the relationships A m A n = e mn k\k must hold, the group multiplication rule 
for parameters xj~ looks 

Xfc — 6rnnkX m X n . (b<2) 

The main claim is that the all properties of any matrix group are straightforwardly deter- 
mined by the bilinear function, the latter is described by structure constants e mn k entering 
the multiplication rule A m A n = e mn kXk- 

• Hamilton-Cayley form [SI [9] [III [111 [T^ rj^ fj^ rj^ [33l [M] • It is possible to expand 
the elements of the group into a power series of linear combination of generators: 

A(a) = i(aiAi H h cln^n), 

because of Hamilton-Cayley theorem this series has three terms for SU (3) and four terms 
for SU(4): 

SU(3), G(a) = e (a)I + ei(a)A(a) + e 2 (a)A 2 (a), 

SU(4), G(a) = e (a)7 + ei(a)A(a) + e 2 (a)A 2 (a) + e 3 (a)A 3 (a). 

2 In the paper we will designate generators in Dirac basis by Ai whereas another set of generators mainly used 
in the literature will be referred as Ai. 
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• Euler-angles representations [ID|l2HI22|l25|l2itl25|l55tH2tl5B|l57tlB5|. In Euler-angles 
representations only a sub-set {A n } C {Xn} of the Lie generators are sufficient to produce 
the whole set (for SU(N) we need only 2 (TV — 1) generators). In that sense all other way 
to obtain the whole set of elements are not minimal. 

In our opinion, we should search the most simplicity in mathematical sense while work- 
ing with basic elements and the structure constants determining the group multiplication 
rule CCO, (USD- 

The material of this paper is arranged as follows. 

In Section [2] an arbitrary 4x4 matrix G G GL(4, C) is decomposed into sixteen Dirac 
matrice^l 



G = AI + iB-/ 5 + iA n l + B a 1 ^ + F mn cr n 



k$ + ka iiq — na 
—Iq — la itiq — ma 



(1.3) 



for definiteness we will use the Weyl spinor basis; four 4-dimensional vectors (k,m,l,n) are 
definite linear combinations of A, B, Ai, Bi, F mn - see (12. 4p . In such parameters (I2.3p . the 
group multiplication law G" = G'G is found in explicit form. 

Then we turn to the following problem: at given G = G(k, m, n, I) one should find parameters 
of the inverse matrix: = G(k' ,mf ,n' ,V) - expressions for [k 1 ,m' ,n' ,V) have been found 

explicitly (for details of calculation see [62J). Also, several equivalent expressions for determinant 
det G have been obtained, which is essential when going to special groups SX(4, C) and its 
subgroups. 

In Section O with the help of the expression for the inverse matrix G^ 1 (k' ,m' ,1' ,n') we 
begin to consider the unitary group SU(A). To this end, one should specify the requirement of 
unitarity G + = G" 1 to the above vector parametrization - so that unitarity conditions are given 
as non-linear cubic algebraic equations for parameters (k, m, I, n) including complex conjugation. 

In Section [4] we have constructed three 2-parametric solutions of the produced equations of 
unitarit}0, these subgroups G\, G2, G3 consist of two commuting Abelian unitary subgroups. 

In Section [5] we have constructed a 4-parametric solutiord - it may be factorized into two 
commuting unitary factors: G = Go (£> SU(2) - see (|5. 15[) . 

The task of complete solving of the unitarity conditions seems to be rather complicated. In 
remaining part of the present paper we describe some relations of the above treatment to other 
considerations of the problem in the literature. We hope that the full general solution of the 
unitary equations obtained can be constructed on the way of combining different techniques 
used in the theory of the unitary group SU (4) and it will be considered elsewhere. 

We turn again to the explicit form of the Dirac basis and note that all 15 matrices are of 
Gell-Mann type: they have a zero-trace, they are Hermitian, besides their squares are unite: 

S P A = 0, (A) 2 = I, (A)+ = A, AG {A fc , k = 1,...,15}. 

Exponential function of any of them equals to 

Uj = e ia i h i = cos aj + i sin ajAj, det & ia ^ = +1, Uf = UJ l , a t G R. 

Evidently, multiplying such 15 elementary unitary matrices (at real parameters Xj) gives again 
an unitary matrix 

JJ — giaiAi e ia 2 A 2 . . . e «ai4Ai4 e ia ;5 Ai 5 jj+ _ e -ia ; A; e ~ia k A k _ _ _ ^-ia^kj ^-iaiki 



3 That Dirac matrices-based approach was widely used in physical context (see [5j [8j |9j [6] 1451 1481 1491 1501 [66] 
and especially [40] V 

4 At this, the unitarity equations may be considered as special eigenvalue problems in 2-dimensional space. 
5 The problem again is reduced to solving of a special eigenvalue problem in 2-dimensional space. 
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At this there arises one special possibility to determine extended Euler angles ai, . . . , ais- For 
the group SU (4) the Euler parametrization of that type was found in [69| . A method to solve 
the problem in [69] was based on the use yet known Euler parametrization for 577(3) - the latter 
problem was solved in [25]. Extension to SU(N) group was done in [70] [TT] . Evident advantage 
of the Euler angles approach is its simplicity, and evident defect consists in the following: we 
do not know any simple group multiplication rule for these angles - even the known solution 
for SU(2) is too complicated and cannot be used effectively in calculation. 

In Section [6] the main question is how in Dirac parametrization one can distinguish SU(3), 
the subgroup in SU(4). In this connection, it should be noted that the basis Aj used in [25] 
substantially differs from the above Dirac basis Aj - this peculiarity is closely connected with 
distinguishing the 577(3) in 577(4). In order to have possibility to compare two approaches we 
need exact connection between Aj and Aj - we have found required formulas^. The separation 
of SL(3, C) in SL(4, C) is given explicitly, at this 3x3 matrix group is described with the help 
of 4 x 4 matrices^- The group law for parameters of SL(3, C) is specified. 

In Section [7] one different way to list 15 generators of GL(4,C) is examined^] 

ai=7°7 2 > a 2 = i7°7 5 , a 3 = 7 5 7 2 , f3 1 =ij 3 j 1 , = i 7 3 , (3 3 = i-f 1 , 

these two set commute with each others ctjflk = Pk&j, and their multiplications provides us with 
9 remaining basis elements of fifteen: 

A 1 =a 1 f3i, B 1 = a 1 (3 2 , C\ = ai/%, 
A 2 = a 2 Pi, B 2 = a 2 (3 2 , C 2 = a 2 f3 3 , 
A 3 = a 3 f3i, B 3 = a 3 f3 2 , C 3 = a 3 f3 3 . 

We turn to the rule of multiplying 15 generators «j, Ai, Bi, Ci and derive its explicit form 
(see CL3D). 

Section [8] adds some facts to a factorized structure of SU(A). To this end, between 9 generators 
we distinguish three sets of commuting ones 

K = {A 1 ,B 2 ,C 3 }, L = {d,A 2 ,B 3 }, M = {B 1 ,C 2 ,A 3 }, 

an arbitrary element from GL(4,C) can be factorized as followg^l 

S = e i3ZjtP e ikK e ilL e irriM^ ^ 

where K, L, M are 3-parametric groups, each of them consists of three Abelian commuting 
unitary subgroups^. On the basis of 15 matrices one can easily see 20 ways to separate SU(2) 
subgroups, which might be used as bigger elementary blocks in constructing a general transfor- 
mational. 

In Sections [9] and [10] we specify our approach for pseudounitary groups SU(2, 2) and 577(3, 1) 
respectively. All generators A' k of these groups can readily be constructed on the basis of the 
known Dirac generators of 577(4) (see ()9.ip ). 



6 This problem evidently is related to the task of distinguishing GL(3, C) in GL{A, C) as well, 
interesting arguments related to this point but in the quaternion approach are given in [35] . 
8 Such a possibility is well-known - see [40j ; our approach looks simpler and more symmetrical because we use 
the Weyl basis for Dirac matrices instead of the standard one as in |40| . 
9 These facts were described in main parts in [40] . 

10 Note that existence of three Abelian commuting unitary subgroups was shown in [40] as well. 
"This possibility was studied partly in [l3l[T4] on the basis of the Hamilton-Cayley approach. 
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2 On parameters of inverse transformations G 1 

Arbitrary 4x4 matrix G G GL(4, C) can be decomposed in terms of 16 Dirac matrices (such an 
approach to the group 7(4, C) was discussed and partly developed in [51 [HI [9l [6l |43 | [4*8 1 [4*9 1 [50l [66] 
and especially in [40]): 



G = AI + iB-f b + iArf + Ba 1 ^ + F mn a mn , 



(2.1) 



where 



7 7+77 



,,ab 



a 



ah 



2g a 

K7V-7V), 



-z 7 °7 1 7 2 7 3 , 



diag(+l, -1,-1,-1). 



Taking 16 coefficients A, B, Ai, B[, F mn as parameters in the group G = G(A,B,Ai,Bi,F mn ) 
one can establish the corresponding multiplication law for these parameters: 



G' = A' I + iflV + iA' a l + B' n l 7 5 + F' mn a mn , 

G = AI + iB-f 5 + iA n l + B n 1 -/ 5 + F mn a mn , 

G" = G'G = A" I + iB" 7 5 + iAU + 7^'V 7 5 + F^ n a mn , 



where 



A" = A' A - B'B - A\A l - B[B l - \F' kl F k \ 
B" = A'B + B'A + A[B l - B[A l + \F^ n F cd e 
A'l = A 1 At - B'Bi + A\A + B[B + A' k F h 



•mncd 



kl 

i_ 171/ a k 1 1 rj/ iTi fcmn , 1 771/ 73 . mrafc 
+ + 2 B k^mn^l + 2 1 mn^k^l , 

B'l = A'Bi + B'Ai - A\B + B[A + B lk F kl 

_i_ J? 1 nfe 1 1 4' 7? c kmn . 1 771/ 4 rank 
^ r lk D + ^^mnt 1 + 2 r mn^ £ Z> 

-^mn = ^ ^im + F mn A ~ {A m A n — A n A m ) — (B m B n — B n B m ) 

1 V R Jkran id* A Jkmn , 1 73/ 771 ^fci , 1 77/ D^fci , / 771/ rpk 771/ rpk 



(2.2) 



The latter formulas are correct in any basis for Dirac matrices. Below we will use mainly Weyl 
spinor basis: 
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a a 
a a 



a a = (I,ai), a a = (I,-an, 
With this choice, let us make 3 + 1-splitting in all the formulas: 

GgGL(4,C), G 



1 



-I 
+7 



ko + ka no — na 
—Iq — la mo — ma 



(2.3) 



where complex 4-vector parameters (k,l,m,n) are defined by [18 



k = A - iB, 
h = Bq — iAq, 



m,Q = A + iB, 
n = B + iA , 



dj + ibj, 



lj = Bj — iAj, no = 7>o + iAg, rij = Bj + iAj. 



(2.4) 



For such parameters (|2,3p . the composition rule (|2.2p will look as follows: 

(k",m";n",l") = (k' ,m';ri ,l')(k,m;n,l), 
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k' ' = k' ko + k'k — ti'qIq + n'l, 

k" = k' k + k'ko + ik' x k — n' l + ti'Iq + in' x Z, 

m'o = m' mo + rn'rn — l' no + Z'n, 

m" = m' m + m'mo — im' x m — l' n + I' no — il' x n. 
n o = k'o n o ~ k'n + riQiriQ + rim, 

n" = k' n — k'riQ + ik' x n + n' m + n'mo — in' x m. 

l'£ = l' Q k + I'k + m' /o - rn'l, 

I" = l' k + I'ko + il' x k + m' l — m'lo — im' x I. 



(2.5) 



Now let us turn to the following problem: with given G = G(k, m, n, I) one should find 
parameters of the inverse matrix: G~ l = G(k' ,m' ,n' ,1'). In other words, starting from 



one should calculate parameters of the inverse matrix G . The problem turns to be rather 
complicated^! the final result is (D = det G, (mn) = mono — mn, and so on) 

k' = D~ 1 [ko(mm) + mo(ln) + lo{nm) — no(lm) + il(m x n)], 
k' = D~ 1 [—k(mm) — m(ln) — l(nm) + n(lm) + 21 x (n x m) 

+ imo(n x I) + il (n x m) + ino{l x m)], 
m'o = -D _1 [fco(/n) + mo(kk) — lo(kn) + n-olk) + m(Z x k)], 
m! = D~ l [-k(ln) - m(kk) + l(kn) - n(kl) + 2nx(lxk) 

+ ino(k x I) + ilo(k x n) + iko(n x I)], 
1' = Z? _1 [+/co(m/) — mo(kl) — lo(km) — no{ll) + im(l x k)], 
I' = D- x [+k(ml) - m(kl) - l(km) - n{ll) + 2mx(kxl) 



G(k, m, n, I) 



+(ko + h) +(h-ik 2 ) 

+(h+ik 2 ) +{k -k 3 ) 

-(lo + k) -ih-ih) 

-{h + ih) ~{h-h) 



+ (n -n 3 ) -(m-OT 2 ) 

-(m + in 2 ) +(n + n 3 ) 

+ (mo-m 3 ) -(mi-im 2 ) 

-(mi + im 2 ) +(m + m 3 ) 



(2.6) 



+ imo(l x k) + iko(l x m) + ilo{m x k)], 
n'o = ko(nm) + mo(kn) — la(nn) — no(km) + ik(m x n)] 

n' = D _1 [-fe(nra) + m(kn) — l(nn) — n(km) + 2k x (m x n) 
+ iko(m x n) + imo(k x n) + ino(m x k)]. 



(2.7) 



Substituting equations (|2.7|) into equation G l G = I one arrives at 



D = k'o = k'oko + k'k — n'oh + w'Z, 

= fc" = fepfc + k'ko + ik' x k- n' l + n'/ + in' x Z, 

I? = m'o = rn'oino + m'm — /^no + Z'n, 

= m" = m' m + m'mo — i m ' x m — l' n + Z'no — iZ' x n. 

= n'o = k'ono — k'n + n' mo + n'm, 

= n" = k' n — k'no + ik' x n + n' m + n'mo — in' x m, 

= I'q = l' Q k + I'k + m' lo - m'l, 

= I" = l' Q k + l'k + H' xk + m' l - m'l - im' x I. 



For more details see [62]; also see a preceding paper [4T] , 
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G 



(2. 



After calculation, one can prove these identities and find the determinant: 

D = det G(k, m, n, I) = (kk)(mm) + (ll)(nn) + 2(mk)(ln) + 2{lk)(nm) - 2(nk)(lm) 
+ 2i[kol(m x n) + mok(n x Z) + lok(n x m) + raoZ(m x k)] 

+ 4(kn)(ml) -4(km)(nl). (2.8) 
Let us specify several more simple subgroups. 

Case A 

Let O-components k$, mo, lo, uq be real-valued, and 3-vectors k, m, I, n be imaginary. Per- 
forming in (|2.5|) the formal change (new vectors are real-valued) 

k =^> ik, rn =^> im, Z iZ, n =^ in, 
ko + ifccr no — ina 

—lo — ila mo — ima 
then the multiplication rules (|2.5p for sixteen real variables look as follows 

k' ' = k k() — k'k — n lo — n'l, 
k" = k' k + k'ko — k' x k — n' l + u'Iq — n' x Z, 
m[j = m' mo — mm — IqUq — I'n, 
m" = m' m + m'rriQ + rn' x m — / n + Z'no + Z' x n, 
n o = ^o n o + k'n + n mo ~~ n'm, 
n" = k' n — k'riQ — k' x n + n' m + n'?no + rt x m, 
Z ' = l' Q k - I'k + m / + m'l, 
Z" = l' Q k + l'k - I' x k + m' l - m'l + m' x Z. 
Correspondingly, expression for determinant (|2.8p becomes 

£> = [tt][mm] + [ll][nn] + 2[m/c][/n] + 2[Ifc][nm] - 2[nk][lm] 

+ 2[/coK m x n ) + mok(n x Z) + l$k{n x m) + riol(m x fe)] 

+ 4(fcn)(mZ) - 4(fcm)(nZ), 

where the notation is used: [ab] = ao&o + 
Case B 

Equations (|2.5p permit the following restrictions: 

and become 

feo = ^o^o + k'k — UqUq + n'n* , 

k" = k' k + k'ko + ik x k — n' n* + n 'tiq + in x n* , 
Uq = k' no — k'n + n' kQ + n'k* , 

n" = k' Q n - k'n + ik' x n + n' fe* + n'/cg - in' x k* ■ 
Determinant D is given by 

D = (kk)(kk)* + (nn)*(nn) + 2(k*k)(n*n) + 2(n*k)(nk*) - 2{nk){nk)* 

+ 2i[kok*(n x n*) — k^k(n* x n) + n* Q n(k x k*) — riQn*(k* x k)] 

+ 4(fcn)(fc*n*) - 4{kk*)(nn*). 
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Case C 

In (|2.9p one can impose additional restrictions 

m o = ko, h = tiq, m = —k. I = —n; (2.10) 

at this G(ko,k,no,n) looks 

_ (k + ika) (no - ina) 
— (no — ina) (ko + ika) 

and the composition rule is 

k' = k' ko — k'k — ti'qTIq + n'n, 
k" = k' k + k'ko — k' x k + n' n + n'no + n' x n, 
n' ' = k' no + k'n + n' ko + n'k, 
n" = k' n — k'riQ — k' x n — n' k + n'ko — n' x k. 

Determinant equals to 

detG = [kk][kk] + [nn][nn] + 2(kk)(nn) + 2(nk)(nk) - 2[nk][nk)) 
+ A(kn)(kn) - 4(kk)(nn). 

Case D 

There exists one other subgroup defined by 
n a = 0, l a = 0, G = 

(mo — ma 

the composition law (|2,5p becomes simpler 

/lq — A^q/cq ~~ I - k k^ k — /lq/u — j - /cq ~\~ %k x 

m o = m o m o + nx'm, rn" = m' m + mlm® — im! x m, 

as well as the determinant D 

detG = (kk)(mm). 

If one additionally imposes two requirements (kk) = +1, (mm) = +1, the Case D describes 
spinor covering for special complex rotation group SO(4, G); this most simple case was consid- 
ered in detail in [18]. 

It should be noted that the above general expression (|2.8|) for determinant can be transformed 
to a shorter form 

detG = (kk)(mm) + (nn)(ll) + 2[kn][ml] 

— 2(kon + nok — ik x n)(mol + lorn + im x I), 

which for the three Cases A, B, C becomes yet simpler: 

(A) : detG = [kk][mm] + [nn][ll] + 2(kn)(ml) 

+ 2(kon + nok + k x n)(mol + fo m — m x I), 

(B) : det G = (kk)(k*k*) + (nn)(n*n*) + 2[&n][fc*n*] 

— 2(kon + nofc — ik x n)(k$n* + n^fe* + ifc* x n*), 

(C) : det G = [A; A;] 2 + [nnf + 2(/cn) 2 - 2(k n + n k + k x nf . 
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3 Unitarity condition 

Now let us turn to consideration of the unitary group SU(4). One should specify the requirement 
of unitarity G + = G^ 1 to the above vector parametrization. Taking into account the formulas 



k% + k*a 



n. 



lq — no m Q 
which can be represented differently 



-l*-l*a 
m*a 



G' 1 



k' + k'a 

-I'n ~ I'd 



n a 



(3.1) 



G + = G{k*,k*;m* ,m*;-%,l*,-n* ,n*), G^ 1 = G(k' ,k';m' ,m';n' ,n' ,l' ,l') 



we arrive at 



k* — k' 
tt — «, , 



7* / 

-t — n , 



k — k . 

I* = n' 



* i/ 

-n — £ , 



rn 
n 



rn 
I'. 



(3.2) 



With the use of expressions for parameters of the inverse matrix with additional restriction 
det G = +1 equations (|3.2p can be rewritten as 



/cq = +ko(mm) + mo(ln) + lo(nm) — no(lm) + il(m x n), 
rriQ = +mo(kk) + ko(nl) + rio(lk) — lo(nk) — in(k x I), 
k* = —k(mm) — m(ln) — l{nm) + n{lm) + 21 x (n x m) 

+ imo(n x I) + i/o(^ x m) + mo(Z x m), 
m* = — m(kk) — k(nl) — n(lk) + l(nk) + 2n x (/ x k) 

— iko{l x n) — ino(l x k) — Uq(ti x fc), 
Iq = +ko(nm) — mo(kn) + lo(nn) + no(km) + ifc(n x m), 

= +rriQ(lk) — fco(mZ) + rio(ll) + lo(mk) — im(l x fc), 
Z* = —k(nm) + m{kn) — l{nn) — n{km) + 2fc x (m x n) 

+ iko(m x n) + im (k x n) + in (m x fc), 
n* = —m(kl) + k{ml) — n(ll) — l(mk) + 2m x (k x I) — imo(k x I) 

— iko(m x I) — il (k x m). 



(3.3) 



Thus, the known form for parameters of the inverse matrix makes possible to write easily 
relations (|3.3|) representing the unitarity condition for group SU(4). Here there are 16 equations 
for 16 variables; evidently, not all of them are independent. 
Let us write down several simpler cases. 

Case A 

With formal changj^l 

k =^> ik, rn =^> im, I il, n in, (3.4) 
equations (|3.3p give 



ko = +ko [mm] + mo[ln] + lo[nm] — no[lm] + l(m x n), 
mo = +mo[kk] + ko[nl] + no[/A;] — lo[nk] — n(k x I), 



3 Let O-components fco, mo, lo, no be real-valued, and 3-vectors k, m, I, n be imaginary. 



12 



V.M. Red'kov, A. A. Bogush and N.G. Tokarevskaya 



k = k[mm] + m[Zn] + l[nm] — n[lm] + 21 x (n x m) 

+ mo(n x I) + Zo(n x m) + no(Z x m), 
m = +m[fc/c] + fe[nZ] + n[lk] — l[nk] + 2n x (I x fc) 

— &o(Z x n) — x k) — lo(n x k), 
Iq = +ko[nm] — mo[kn] + lo[nn] + no[&m] + k(n x m), 
riQ = +rrio[lk] — ko[ml] + no[ll] + Zo[mA] — m(l x k), 
I = +k[nm] — m[kn] + l[nn] + n[km] + 2k x (m x n) 

+ fco(m x n) + mo(fe x n) + no(m x fe), 
n = +m[kl] — k[ml] + n[ll] + l[mk] + 2m x (k x I) 

— mo(k x I) — ko(m x I) — l (k x m). 

Here there are 16 equations for 16 real-valued variables. 

Case B 

Let 

mo = fcg, m = k* , Iq = rig, I = n* , 

ko = TOg, fe = m* , no = Iq, n = I* , 

or symbolically m = k* , I = n* . The unitarity relations become 

kl = +k (k*k*) + kl(n*n) + n* Q {nk*) - n (n*k*) + in*(k* x n), 
k* = -k(k*k*) - k*(n*n) - n*{nk*) + n(n*k*) 

+ 2n* x (n x k*) + ik^n x n*) + in* Q (n x k*) + in (l x m), 
Uq = +kl(n*k) - k (k*n*) + n (n*n*) + nl(k*k) - ik*(n* x k), 
n* = -k*(kn*) + k{k*n*) - n(n*n*) - n*{k*k) 

+ 2k* x (k x n*) — ik^ik x n*) — iko(k* x n*) — in^k x k*), 

and 8 conjugated ones 

ko = +k*i(kk) + ko(nn*) + no(n*k) — n* Q {nk) — in{k x n*), 
k = -k*(kk) - k(nn*) - n(n*k) + n*(nk) 

+ 2n x (n* x k) — iko{n* x n) — ino{n* x k) — in^n x k), 
no = +ko(nk*) — ko(kn) + n^(nn) + no{kk*) + ik(n x k*), 
n = —k(nk*) + k*(kn) — n*{nn) — n(kk*) 

+ 2k x (k* x n) + iko(k* x n) + iko(k x n) + ino(k* x k). 

It may be noted that latter relations are greatly simplified when n = 0, or when k = 0. Firstly, 
let us consider the case n = 0: 

k* = +k (k*k*), k* = -k(k*k*). 
Taking in mind the identity 

detG= (kk)(kk)* = +1 =^ (kk) = +l, (kk)* = +l, 
we arrive at feg = +ko, k* = —k. It has sense to introduce the real-valued vector c a : 

&0 = +&o = Co, k* = —k : k = ic, 
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then matrix G is 

G(k,m = k*, 0,0) = c ° + tca °_ ~SU(2). 

Co - ico 

Another possibility is realized when k = 0: 

Uq = +no(nn)*, n* = —n(nn)*. 
With the use of identity 

detG = (nn)(nn)* = +1 =^ (nn) = +1, (nn)* = +1, 
we get 

n* Q = +no = Co, n* = — n, n = ic, 
Corresponding matrices G(0, 0,n, I = n*) make up a special set of unitary matrices 



G 



co — ico 

(cq + ico) 



G 



— (co — ico) 

(cq + ico) 



(3.5) 



However, it must be noted that these matrices ()3.5[) do not provide us with any subgroup because 
G 2 = -I. 



Case C 

Now in equations (|3.4|) one should take 
mo = &o, Zo = n o 5 m = —k, 

then 



-n, 



ko = +k$[kk\ + A;o(nn) + rio(nk) — rio[nk], 

k = k[kk] — k(nn) — n(nk) — n[nk] + 2n x (n x fc), 

no = +ko(nk) — k^lkn} + no[nn] + rio(kk), 

n = —k(kn) — k[kn] + n[nn] — n(kk) + 2k x (fc x n). 



(3.6) 



4 2-parametric subgroups in SU (4) 

To be certain in correctness of the produced equations of unitarity, one should try to solve them 
at least in several most simple particular cases. For instance, let us turn to the Case C and 
specify equations (|3.6|) for a subgroup arising when k = (ko, fci,0,0) and n = {nQ,n\,0, 0): 

ko = +ko[kk] + ko(nn) + no(nk) — rio[nk], 
k\ = +k\[kk] — k\(nn) — n\(nk) — n\[nk], 
no = +ko(nk) — ko[kn] + no[nn] + rio(kk), 

rii = —ki(kn) — k\[kn] + ni[nn] — rii(kk), (4.1) 



they are four non-linear equations for four real variables. It may be noted that equations (|4.1 
can be regarded as two eigenvalue problems in two dimensional space (with eigenvalue +1): 



^ + nl)-l + {kl-n\) 
— 2n\k\ 



'4 + nl) 



-2n\k\ 

-l-{kl-nl) 





k Q 









no 
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l + n{)-\ + {kl-nl) 


-2n k 




ki 







-2n k 


{k\ + nf) - 1 - {kl - n§) 




ni 








The determinants in both problems must be equated to zero: 

[{kl + nl)-lf-{kf-nif 
[(kf + nl) - if - {k 2 - n 2 ) 



An\k\ 
4nX 



0, 
0, 



or 



0. 



+ n 2 ) - I] 2 - {kl + nl? = 0, [{kl + n\) - if - {k 2 + n 2 ) 2 

The latter equations may be rewritten in factorized form: 

[{kl + n 2 ) - 1 - {kl + nl)][{k 2 + nl) - 1 + {kl + n 2 )] = 0, 
[{kl +nl)-l- {kl + n 2 )][{kl +nl)-l + (A; 2 + n 2 )] = 0. 

They have the structure: AC = 0, BC = 0. Four different cases arise. 
(1) Let C = 0, then 



kl + nl + kl + nl = +1. 

(2) Now, let A = 0, 5 = 0, but a contradiction arises: A + B = 0, A + 5 

(3) — (4) There are two simples cases: 



-2. 



A = 0,C 

5 = o,c 








kl + n-Q 
kl + ?i 2 



fcl = 0, ni 
A; = 0, n 



0, 
0. 



Evidently, (|4.3p and (|4,4p can be regarded as particular cases of the above variant (|4,2 
one should take into account additional relation 



1, 



detG = [kk][kk] + [nn][nn] + 2{kk){nn) 

+ 2{nk){nk) - 2[nk][nk]) + A{kn){kn) - A{kk){nn) 

which can be transformed to 

det G ={kl + kl + nl + nl) 2 - A{hn + k m) 2 = +1. 
Both equations (|4.2p and (|4.5p are to be satisfied 

{kl + n 2 Q + kl + nl) = 1, {kl + kl + + n 2 ) 2 - 1 - A{k x n + k nif = 0, 
from where it follows 

k\UQ + k n\ =0, kl +nl + kl+nl = +1. 

They specify a 2-parametric unitary subgroup in 5 £7 (4) 

Gf = G7 1 , detGi = +l, 

/c 2 + + /c 2 + nl = +1, 



fcin + k ni = 0, 

fco + ikicr 1 no — inia 1 



Gx 



-no + in\/j l ko + ik\a l 



ko 


ik\ 


n 


—in\ 


ik\ 


ko 


—in-y 


no 


-n 


ini 


ko 


iki 


in\ 


-n 


iki 


ko 



(4.2) 



(4.3) 
(4.4) 

Now, 



(4.5) 



(4.6) 
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Two analogous subgroups are possible: 

G+ = G^\ detG 2 = +l, 

k 2 n + k n 2 = 0, k% + n^ + k\ + n 2 = +1, 



Go 



feo + ik 2 o 2 tlq — iri2(J 2 
-no + in 2 a 2 ko + ife 2 cr 2 





fc 2 


»o 


-n 2 


-k 2 


ko 


n 2 


»o 


-n 




ko 


fc 2 


-n 2 


-n 


-k 2 


ko 



(4.7) 



G^Gg 1 , detG 3 = +l, 

fe 3 n + k n 3 = 0, feg + n + ^3 + n 3 = + 1 ' 



G, 



(fe + ife 3 ) 

(k - ik 3 ) 

-(n -in 3 ) 

-(n + in 3 ) 



(n - in 3 ) 


(k + zfc 3 ) 






(n + m 3 ) 


(k - ik 3 ) 



(U 



Let us consider the latter subgroup 
ters is 

fed' = k k - k' 3 k 3 - n' n + n' 3 n 3 , 
n'o = fed n o + k' 3 n 3 + n k + n' 3 k 3 , 



in some detail. The multiplication law for parame- 

k 3 = k' k 3 + k 3 k + n' n 3 + n' 3 no, 
n 3 = k' Q n 3 - k' 3 n - n' Q k 3 + n' 3 k . 



For two particular cases (see f)4.3[t and (|4.4p ). these formulas take the form: 



{G } 



fe| + n| = l, 


feo 


= 0, 


n = 0, 


k'o = ~k' 3 k 3 +n' 3 n 3 , 


^'3 


= 0, 




n = +k' 3 n 3 + n' 3 k 3 , 




= 0, 




^0 + n = 1' 


k 3 


= 0, 


n 3 = 0, 


fed' = fed^o - n' no, 








n'o = fed n o + ^d^o- 









(4.9) 



Therefore, multiplying of any two elements from G 3 ° does not lead us to any element from 
G 3 °, instead belonging to G°: G 3 G 3 G G°. Similar result would be achieved for G\ and G 2 : 
G\ G\ G G° , G2 G° G G° . In the subgroup given by (|4.9p one can easily see the structure of the 
1-parametric Abelian subgroup: 



feo = cos a, 
G°(a) = 



n 


= sin a, 






cos a 





sin a 








cos a 





sin a 


— sin a 





cos a 








— sin a 





cos a 



a G [0,27r]. 



(4.10) 



In the same manner, similar curvilinear parametrization can be readily produced for 2- 
parametric groups (|4.6p - (|4.8p . For definiteness, for the subgroup G 3 such coordinates are 
given by 

feo = cos a cos p, k 3 = cos a sin p, 

no = sin a cos p, — n 3 = sin a sin p, a G [0, 2n] , 
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and matrix G-? is 



G 3 (p,a) 



cos ae 


— sin ae ip 




sm ae 



cos ae 




-i>p 



smae 



-ip 



cos ae 




ip 



smae 




cos ae 



->i> 



-ip 



(4.11) 



One may note that equation (|4.1ip at p = will coincide with G°(a) in (|4,10p : G 3 (p = 
0, a) = G°(a). Similar curvilinear parametrization may be introduced for two other subgroups, 
G\ and G 2 . 

One could try to obtain more general result just changing real valued curvilinear coordinates 
on complex. However it is easily verified that it is not the case: through that change though 
there arise subgroups but they are not unitary. Indeed, let the matrix (14. 10h be complex: then 
unitarity condition gives 



cos a cos a + sm a sm a 



cos a sm a + sm a cos a 



0. 



These two equations can be satisfied only by a real valued a. In the same manner, the the 
formal change {G±,G2,Gs} ==> {G*p, G^, G^} again provides us with non-unitary subgroups. 

It should be noted that each of three 2-parametric subgroup G\, G2, G 3 , in addition to Go (a), 
contains one additional Abelian unitary subgroup: 



K 3 



ko + ik\a l 

ko + ik\a l 

ko + i/c 2 <7 2 

k + ikio 2 



ko + ik 3 a 3 







k + ik 3 cr 3 



Kx C Gi, 
K 2 C G 2 , 
K 3 C G 3 , 



k^ + k\ — 1, 



^0 + ^2 — 1) 



1. 



It may be easily verified that 

Gi = Go-f^i = K\Gq, G2 = GqK 2 = K 2 Gq, G3 = GqK 3 = K 3 Gq. 
Indeed 

cos ako + i cos akia 1 sin afco + i sin ak\ 



Go{a)Kx = KxGoia) 



sin ako — i sin ak\ cos afco + i cos ak\o 



and with notation 

cosafco = k'o, cosaki = k[, sinafco = n , 
k'^n^ + k'xn'o = 0, k 2 + k 3 + n 2 + n 3 2 = 1 

we arrive at 

k' Q + ik^a 1 n' — in\ 
—n' + k' + z'A^er 1 



sin afci 



-n 



G0K1 = KiGq 



C Gi. 



5 4-parametric unitary subgroup 

Let us turn again to the subgroup in GL(4, G) given by Case C (see (|2.10p ): 



G 



(ko + k&) (no — na) 
(no — na) (ko + ko) 
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when the unitarity equations look as follows: 



ko = +ko[kk] + ko(nn) + no(nk) — no[nk], 

riQ = +ko(nk) — k^[kn] + no[nn] + no(kk), 

k = k[kk] — k(nn) — n(nk) — n[nk] + 2n x (n x k), 

n = —k(kn) — k[kn] + n[nn] — n{kk) + 2k x {k x n). 



They can be rewritten as four eigenvalue problems: 

(+1) 



ko 
n 



[kk] + (nn) ink) — [nk] 
(nk) — [nk] (kk) + [nn] 

+ ([kk] - [nn]) -2{nk) 
-2(nk) -{[kk] - [nn]) 

+ ([kk] - [nn]) -2{nk) 
-2{nk) ~{[kk] - [nn]) 

+ ([kk] - [nn]) -2(nk) 
-2(nk) ~([kk] - [nn]) 

These equations have the same structure 



k 
no 



(5.1) 





ki 


= (+1) 


h 




ni 


ni 




k 2 


= (+1) 


k 2 










h 


= (+1) 


k 3 




n3 


^3 



(5.2) 



A 


C 




Zi 


= A 


Zi 


C 


B 




z 2 


z 2 



where A = +1. Non-trivial solutions may exist only if 



det 



A- A C 
C B- A 



which gives two different eigenvalues 



A I 



A + B + y/(A - B) 2 + AC 2 



A + B + y/{A- B) 2 + 4C 2 



In explicit form, equations (|5.1|) looks as follows: 
A 



A C 
C B 



ko 
no 



k 
n 



A 



k z + ng) + (k z - n z ), B = (k z + ng) - (k z - n% 
\y = (k 2 + n 2 ) + J(k 2 -n 2 ) 2 + 4(kn) 2 , 



C 



2kn, 



(5.3) 



A 2 = (k 2 + nl) ~ ^(k 2 -n 2 ) 2 + A(kn) 2 . 
The eigenvalue A = +1 might be constructed by two ways: 



Ai=+1, k 2 +n 2 = l-yJ(k 2 -n 2 ) 2 +A(kn) 2 , 
A 2 = +l, k 2 + n 2 = 1 + ^(k 2 - n 2 ) 2 + A(kn) 2 . 
These two relations (|5,4h are equivalent to the following one: 
(l-k 2 -n 2 ) 2 = (k 2 -n 2 ) 2 + A(kn) 2 . 



(5.4) 
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Thus, equations (|5.3p have two different types: 
Type I 

(A - l)ko + Cn = 0, Ck + (B — l)n = 0, 

k 2 Q + nl = 1 - ^(k 2 - n 2 ) 2 + 4(fcn) 2 , (5.5) 
k% + n 2 <+l, {k 2 - n 2 ) 2 + A{kn) 2 < +1; 

Type II 

(A - l)k + Cn = 0, Ck + {B - l)n = 0, (5.6) 
fc 2 + ng = 1 + ^(fc 2 -n 2 ) 2 + 4(fcn) 2 , A; 2 , + n 2 > +1. 

Now let us turn to equations (|5.2p . They have the form 



A 


C 




ki 


= A 


ki 


C 


-A 




m 


m 



i = 1,2,3, 



A = fcg + k 2 - nl - n 2 , B = -A, C = -2(k n - kn), 

Ai = +\J {kl + k 2 -n 2 - n 2 ) 2 + A{k n - kn) 2 , 

A 2 = - ^{k 2 + k 2 - n 2 - n 2 ) 2 + 4(A; n - kn) 2 . 

As we are interested only in positive eigenvalue A = +1, we must use only one possibility 
A = +1 = Ai, so that 

(A-l)k + Cn = 0, Ck - (A + l)n = 0, (5.7) 
1 = [kl + k 2 - nl - n 2 ) 2 + A(k n - kn) 2 . 

Vector condition in (|5.7p says that k and n are (anti)collinear: 

k = Ke, n = Ne, e 2 = 1, e 6 5 2 , (5.8) 

so that (|5.7p give 

(A - l)liC + CiV = 0, - (A + 1)7V = 0, 

1 = (fcg + K 2 - n 2 - N 2 ) 2 + 4(fc n - KN) 2 , (5.9) 
A = kl + K 2 - nl - N 2 ,C = -2(k n - KN). 

With notation (|5.8p . equations (|5.5p - (|5.6p take the form: 
Type I 

(A - l)k + Cn = 0, C^o + (5 - l)n = 0, fc 2 + ?ig = 1 - (K 2 + A^ 2 ), 
Type II 

(A - l)k + Cn = 0, C&o + (5 - l)no = 0, k 2 Q + n 2 = I + {K 2 + N 2 ), (5.10) 

where 

A = (kl + n 2 ) + (K 2 -N 2 ), B={kl + nl)-{K 2 -N 2 ), C = —2KN. (5.11) 
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Therefore, we have 8 variables e, ko, no, K, N and the set of equations, (|5.9p - (|5.1ip for them. 
Its solving turns to be rather involving, so let us formulate only the final result: 



h, 

'.2 

G 



k = Ke, 



k 2 + K 2 + n 2 + N 2 



no, n = Ne, 

-riy -+1, k N + n K = 0, 

(h + iKea) (no — iNea) 
(no — iNea) (ko + iKea) 



(5.12) 



~ 2 + K 2 + nl + N 2 ) 2 = +1. 



It should be noted that 

det G — Tn t "-o 
The unitarity of the matrices (|5.12p may be verified by direct calculation. Indeed, 

G + 



(h — iKea) — (no + iNea) 
(no + iNea) (ko — iKea) 



and further for GG + = I we get (by 2x2 blocks) 

(GG + )u = k 2 + K 2 + n 2 + N 2 = +1, (GG + ) 12 = -2i(n K + A; iV)(e«T) = 0, 
(GG + ) 22 = fcg + K 2 + nl + N 2 = +1, (GG + ) 2 i = +2i(n K + A; iV)(e«T) = 0. 

One different way to parameterize (|5.12|) can be proposed. Indeed, relations (|5.12|) are 



ko, 



k = Ke, 



no, 



n = Ne, 



K 2 N 2 
fc o(l + T^)+^o(l + —) = +!, 

K n TLr. 



K 

h 



N 
no 



W, 



or 



n = —noWe, 
K = k W, N 



ko, k = hWe, no, 
(k 2 +n 2 )(l + W 2 ) = +l, 

Therefore, matrix G can be presented as follows: 

k (l + iWea) n (l+iWea 

-no(l + iWea) ko(l + iWea 

+ n 2 )(l + W 2 ) = +1 



-n Q W, 



< kl + n\ < 1. 



G 



(5.13) 



W 



k 2 + nl 



Evidently, it suffices to take positive values for W. The constructed subgroup (|5.13p depends 
upon four parameters ko, no, e: 

0<fcg + ng<l, e 2 = l, (k 2 + n 2 )(l + W 2 ) = +1. 

Let us establish the law of multiplication for four parameters ko, no, W = We: 

G" = 



k' (l + iW'a) n' (l + iW'a) 
-ri (l + iW'a) k' (l + iW'a) 



k (l + iWa) n (l + iWa) 
-n (l + iWa) k (l + iWa) 



or by 2 x 2 blocks 

(11) = (k' k - n n )(l + iW'a)(l + iWa), 
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(12) = (k' n + n' k )(l + iW'd)(l + iWa), 

(21) = -{k' n + n' k )(l + iW'a)(l + iWa), 

(22) = (k' k - n' n )(l + iWa){\ + iWa). 

As (11) = (22), (12) = —(21); further one can consider only two blocks: 

/ w' + W - W' x W 

(11) = (k' k - n' n )(l - WW) \l + i + _ w , w 3 

/ W' + w _ w' x w ' 

(12) = (£> + n' k )(l - W'W) 1 + i — ...... a 



1 - WW 



So the composition rules should be 



ftg = (k' k - n' n )(l - W'W), n' ' = (k' n + n' fc )(l - W'W), 
W + W -W X W 



w" 



1- WW 



The later formula coincides with the Gibbs multiplication rule (see in [32]) for 3-dimensional 
rotation group SO(3,R). It remains to prove the identity: 

(fc 2 + n 2 )(l + W 2 ) = +l 

which reduces to 

(k' k - n' n ) 2 + (k' n + n^o) 2 [(1 - W'W) 2 + (W' + W - W x W) 2 ] . (5.14) 
First terms are 

{k' k - n' n ) 2 + (k' n + n' k ) 2 = (k' 2 + n£)(kl + n 2 ,). 
Second term is 

(1 - W'W) 2 + (W + W - W' x W) 2 = (1 + W 2 )(l + W 2 ). 
Therefore, (I5.14|) takes the form 

(k'o 2 + n' 2 ){k 2 + n 2 )(l + W' 2 )(l + W 2 ) = l 
which is identity due to equalities 

(k' 2 + n' 2 )(l + W' 2 ) = 1, (k 2 + n 2 )(l + W 2 ) = 1. 

It is matter of simple calculation to introduce curvilinear parameters for such an unitary 
subgroup: 



e = (sin cos 0, sin sin ^, cos 0), 

ko = cos a cos p, K = cos a sin p, no = sin a cos p, 



N 



sin a sin p, 



and G looks as follows 



A £ 
-S A 



cos a(cos p + i sin/? cos 0) i cos a sin p sin 0e 11 
i cos a sin p sin Qe 1 ^ cos a (cos p — i sin p cos I 

sin a(cos p + i sin p cos 0) +i sin a sin p sin 9e~ l< ^ 
i sin a sin p sin 9e 1 ^ sin a(cos p — i sin p cos 0) 
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It should be noted that one one may factorize 4-parametric element into two unitary factors, 
1-parametric and 3-parametric. Indeed, let us consider the product of commuting unitary groups, 
isomorphic to Abelian group Go and SU(2): 



G = G ® SU(2) = SU(2) <g> G 



k' 

-n'n 



V 



k' ao + ik' Q aa 



-n Q a 
with the notation 



n' ao + in' aa 
k' Q ao + ik' aa 



ao + iaa 




k' 2 4-T7 /2 





ao + iao 



1, 



ag + a' 



1, 



fc a = &o, 

( fc 2 +n 2 )(l + VF 2 



k W, n' a = n , 
2 ag + n 2 ag)(l + W 2 ) 



J2. 
"0 



n' n a 



(A- 



n' 2 )(a 2 + a 2 ) 



takes the form 



G <8> 5C/(2) = 5C/(2) G 



^(l+i^ff) n (l + iWa) 
-noJl + iWa) /c (l + iW<r) 



Let us summarize the main results of the previous sections: 



G. 



(5.15) 



Parametrization of 4 x 4 matrices G of the complex linear group GL(4, G) in 
terms of four complex vector-parameters G = G(k, m, n, I) is developed and the 
problem of inverting any 4x4 matrix G is solved. Expression for determinant of any 
matrix G is found: det G = F(k, m, n, I). Unitarity conditions have been formulated 
in the form of non-linear cubic algebraic equations including complex conjugation. 
Several simplest solutions of these unitarity equations have been found: three 2- 
parametric subgroups Gi, G2, G3 - each of subgroups consists of two commuting 
Abelian unitary groups; 4-parametric unitary subgroup consisting of a product of 
a 3-parametric group isomorphic SU(2) and 1-parametric Abelian group. 

The task of full solving of the unitarity conditions seems to be rather complicated and it will 
be considered elsewhere. In the remaining part of the present paper we describe some relations 
of the above treatment to other considerations of the problem in the literature. The relations 
described give grounds to hope that the full general solution of the unitary equations obtained 
can be constructed on the way of combining different techniques used in the theory of the unitary 
group 577(4). 



6 On subgroups GL(3, C) and SU(3), expressions 
for Gell-Mann matrices through the Dirac basis 

In this section the main question is how in the Dirac parametrization one can distinguish 
GL(3,C), subgroup in GL(4,C). To this end, let us turn to the explicit form of the Dirac 
basis (the Weyl representation is used; at some elements the imaginary unit i is added) 



17 



-10 
0-100 

10 

1 

-i 

-i 

i 

1 
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7 y 



0010 
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0010 
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2(7 



01 



2ia 



12 



-i 

i 

-i 

1 

10 
10 
-1 
0-10 

10 
0-100 
10 
-1 



2a 



02 



2ia 



23 



-i 

i 

1 
-i 

-t 

1 
i 
-i 

10 
10 
1 
10 



7 5 7 3 



2(7 



2ia 



31 



10 
-1 
10 
0-100 

10 
0-100 
0-10 


-i 

i 








1 

-i 



-% 

1 



All these 15 matrices Aj are of Gell-Mann type: they have a zero-trace, they are Hermitian, 
besides, their squares are unite: 



A, 



SpA = 0, (A) 2 = I, (A)' 
Exponential function of any of them equals to 



A€{A fc :*; = !,..., 15}. 



U 



JaA 



cos a + i sin a A, 



det e 



iaA 



+ 1, 



u 



a £ R. 



Evidently, multiplying of such 15 elementary unitary matrices (at real parameters Xj) results in 
an unitary matrix 



U 



e ia 1 A le ia 2 A2 



e iai 4 Ai4 e ia i5 A( 5 



At this there arise 15 generalized angle- variables a±, . . . ,015- Evident advantage of this ap- 
proach is its simplicity, and evident defect consists in the following: we do not know any simple 
group multiplication rule for these angles. 

It should be noted that the basis Aj used in [69] substantially differs from the above Dirac 
basis Aj - this peculiarity is closely connected with distinguishing SU (3) in SU(4). This problem 
is evidently related to the task of distinguishing GL(3, C) in GL(4, C) as well. 

In order to have possibility to compare two approaches we need exact connection between Aj 
and Aj. In [69] the following Gell-Mann basis for SU(4) were used: 



vio 
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1 

10 



U4 







-i 

i 



\15 



1 



10 
10 
10 
-3 



All the A excluding As, A15 possess the same property: 
A? = +Ai, i + 8, 15. 
The minimal polynomials for As, A15 can be easily found. Indeed, 



10 

10 

4 





(As) 3 





1 











1 





1 








3^ 








-8 




















therefore 

(As) 3 = § As - -^(As) 2 - 
Analogously, for A15 we have 



(A 



15 J 



10 
10 
10 
9 





1 











1 





1 
















1 
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and 



^Ai5 



^(Ais) 



(A15) - 2 ~ i0 ^ 
Comparing Aj and A, one can readily derive the linear combinations: 
7 ° + 7 5 7 3 = 2A 4 , 7°-7 5 7 3 = 2An, i 7 5 7 ° + »7 3 = 2A 5 , 



i 7 5 7 ° 



^ 7 



i 7 i + 7 5 7 2 
,01 



23 



" 2Ai2, 

2Aio, 
= -2A13, 



7V + «7 2 



1 



Z7 



2cr ul - 2ia 
and additional six combinations 



7 V - iT 2 = 2A 9 , 
2a m + 2ia 23 = 2Ai, 
2cr 02 + 2io- 31 = 2A 2 , 2ct 02 - 2za 31 = 



7 5 7 2 



- 2A6, 
2A 7 , 



-2A 



14, 



2o- 03 + 2ia 



7 5 + 2a 03 



7 5 + 2ia 12 



12 





-2 







2a 



03 



2ia 



12 


















-2 





























2 

















-2 














2 


















2a 1 



03 



2io 



12 







0-2 



2 



2 



-2 





2 







+2 



(6.2) 
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they should contain three linearly independent matrices. Those three linearly independent 
matrices might be chosen in different ways. Let us introduce the notation: 



A 



10 
0-100 







0-1 

+1 



B 





-10 

1 



-10 





1 



c 



-10 



10 




0-100 



1 



The matrices a, b, c have the 3x3 blocks different from zero, so they could be generators 
for SU(3) transformations; whereas A, B, C may be generators only of the group SU(4:). All 
six matrices a, b, c, A, B, C have the same minimal polynomial: 

A 3 = A. 

Linear space to which these six matrices a, b, c, A, B, C belong is 3-dimensional. Indeed, one 
easily obtains 



c = b-a, C - A 
basis {a, b, C} = 



b, C - B = a, B-A- 
c = b — a, A = C — 6, 



B = C-a. 



These relations can be rewritten differently 



a = b- c, C - A 
basis {b, c, A} == 



b, C-B = b-c, B-A = c, 

a = b- c, C = A + b, B = A + c 



or 



b = a + c, C - A 
basis {a, c, B} = 



b, C-B = b-c, B-A = c, 
b = a + c, C = B + a, A = B 



One should note that in the basis Aj (|6.2p the corresponding three linearly independent 
elements A3, As, A15 are taken as follows: 



A 3 = 
A15 



i2ff 03 + i2i<r 12 , 



^-^ 5 + 3^ ra -37J 2! " 12 - 



1 „,5 1 n„03 1 1 r^ 12 

Ve 7 "75 2cJ + VE 2ia ' 



their minimal polynomials look 



(A3) 3 = A 3 , 



(As) s 



rAv 



(A 



15- 



^Ai 5 



It is the matter of simple calculation to find relationships between A3, As, A15 and the basis 
{a,b,C} : 



A3 — a, As — -^a — -^b, A15 - 
In the following we will use the notation 
a = A3, b = A'g, C = A' 15 , 



(6.3) 
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so the previous formulas (|6.3p will read ({A3, As, A15} -4=>- {A3, A' 8 , A' 15 }) 
A 3 = A 3 , A 8 = -±X 3 - -J=A' 8 , A15 = ^A 3 + -j=A 8 - ^gA' 15 . 



The inverse relations are 

A3 = A3, A' 8 = jX 3 - ^A 8 , A' 15 = — 3^Ai5 + |A3 — 2-^Ag. (6.4) 

Now, starting with the linear decomposition of G E GL(4, C) the in Dirac basis (|2.ip : 

G = a I + ib j 5 + iA j° + iA kl k + £ 7°7 5 + #fc7 V 
+ a k 2a 0k + h2a 23 + 6220-31 + & 3 2<7i 2 
= a I + ib j 5 + iA j° + A k {i-/ k ) + iB (i~f 5 j°) - B k (j 5 j k ) 
+ a k {2aQ k ) - ibi(2ia 23 ) - ib 2 (2ia 3l ) - ib 3 (2ia u ), 

with the help of the formulas 

7 + 7 5 7 3 = 2A4 1 7 _ 7 5 7 3 = 2An ^ - 7 5 7 + - 7 3 = 2 ^ > 

i 7 5 7 ° - = 2A12, 7V + i7 2 = 2A 6 , 7V - i 7 2 = 2A 9 , 

i 7 i + 7 5 7 2 = 2A10, h 1 ~ 7 5 7 2 = 2A 7 , 2cr 01 + 2ia 23 = 2Ai, 

2cr 01 - 2ia 23 = -2Ai 3 , 2cr 02 + 2ia 31 = 2A 2 , 2o- 02 - 2ia 31 = -2\ u , 

2a m + 2ia 12 = 2A 3 , 2a 03 - 2ia 12 = 2(A' 15 - A' 8 ), 7 5 + 2ct 03 = 2A' 15 , 

7 5 -2a 03 = 2(A 8 -A 3 ), 7 5 + 2^ 12 = 2A' 8 , 7 5 - 2ia 12 = 2(A' 15 - A3) 

and inverse ones 

7° = A 4 + Aii, 7 5 7 3 = A 4 -An, i7 5 7° = A5 + A12, = A 5 - A i2 , 

7 5 7 X = A 6 + A 9 , «7 2 = A 6 - A 9 , fy 1 = X w + A 7 , 7 5 7 2 = Ai - A 7 , 
2k 23 = Ai + A 13 , 2a 01 = Ai - A13, 2icj 31 = A2 + A14, 2a 02 = A 2 - A i4 , 
2a 03 = A 3 - (A' 8 - A' 15 ), 2m 12 = A 3 + (A' 8 - A' 15 ), 7 5 = -A3 + (A' 8 + A' 15 ), 

we will arrive at 

G = a I + (ai - i6i)Ai + (a 2 - ib 2 )X 2 + {iA - B 3 )X 4 + (A 3 + iB )X 5 + {A 2 - £i)A 6 
+ (Ai + B 2 )X 7 + (a 3 - ib 3 - ib )X 3 + (-ib 3 + ib - a 3 )X' 8 + (ib + a 3 + i6 3 )A' 15 
+ (-£1 - A 2 )A 9 + (Ai - B 2 )X 10 + (iA + 5 3 )An + (-A 3 + iB )X 12 
+ (-ai - ifei)Ai 3 + (-a 2 - ib 2 )X u . 

In variables (k,m,l,n) (see ((231), flZ3])) 

Bq — iAq = Iq, Bj — iAj = lj, Bo + iAq = no, Bj + iAj = rij, 
Oo — ibo = ko, aj — ibj = kj, ciq + ibo = uiq, cij + ibj = rrij 

the previous expansion looks 

G = i(&o + m )I + fciAi + /c 2 A 2 + |[(n - n 3 ) - (Z + Z 3 )]A 4 + ^[-(n - n 3 ) - (Z + Z 3 )]A 5 
+ ^[-( n i + in i) ~ (h ~ ih)]Xe + ji[{ni + in 2 ) - (h - il 2 )]X 7 + [k 3 + ±(/c - ^0)] A3 
+ [~m 3 + i(m - k )] A' 8 + [m 3 + \{m - k )] X' 15 + |[-(ni - m 2 ) - (Zi + iZ 2 )]A 9 
+ ^K n i ~ in 2) - ('1 + ^2)]Aio + \[{no + n 3 ) - (l - l 3 )]Xn 
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+ 4[-(«o + "-3) - (lo ~ k)]^i2 ~ m 1 X 13 - m 2 Ai 4 . 



Let coefficients at A9, A10, An, A12, A13, A14 be equal to zero: 

(ni - in 2 ) + {h + ih) = 0, (ni - m 2 ) - (h + ih) = 0, 
(n + n 3 ) + (Iq - h) = 0, m 1 = 0, m 2 = 0. 

Note that we do not require vanishing of the coefficient at A' 15 : 
A i5 m 3 + = + 0. 



(n + n 3 ) - (l - h) = 0, 
(6.5) 



As a result we have a subgroup of 4 x 4 matrices defined by 10 complex parameters. At this 
four elements are diagonal matrices: 



Aa 



A' 



I, A3, /Vg, /A 15 , 

all other matrices have on the diagonal only zeros. Equations (|6.5p give 

m 2 = ni, n 3 = -n , il 2 = —h, h = mi = 0, m 2 = 0, 

so that any matrix G(k a , n , ni, Zq^i, m , m 3 ) is decomposed according to 



(6.6) 



G = k\\i+ k 2 X 2 + (n - /o)A4+ i(no + h)^5 + (— n\— h)X& + i(-n\+ h)\ 7 + \(k Q + mo) I 
+ [k 3 + i(/eo - m )] A 3 + [-m 3 + ^(m - feo)] A 8 + ["i3 + - ko)] A' 15 . (6.7) 



Explicit form of the matrices parameterized by (|6.7p can be obtained from representation for 
arbitrary element of GL(4, C) (|2.6|) 



G(fc, m, n, I) 



+ (k + k 3 ) +(k 1 -ik 2 ) + (n -n 3 ) -(m-in 2 ) 
+ (/ci+ifc 2 ) +(A; -fe 3 ) -(m+m 2 ) +(n + n 3 ) 
il 2 ) +(m -m 3 ) -(mi-im 2 ) 



-(lo + h) 
-(h+ih) 



ih 

-(lo-h) ~(m 1 + im 2 ) + (m + m 3 ) 



with additional restrictions (16.61): 



G 



ko + k 3 k\ — ik 2 +2no 

k\ + ik 2 ko — k 3 —2n\ 

— 2/o — 2Zi mo — m 3 

mo + m 3 



If additionally one requires mo + m 3 = 1, then 



ko + k 3 k\ — ik 2 +2no 

k\ + ik 2 ko — k 3 —1rt\ 

-2/ -2/i 1 - 2m 3 

1 



(6.8) 



G 



with decomposition rule 

G = k\\\ + k 2 \ 2 + (no - Zo)A4 + i(n + h)^5 + (—ni - h)^6 + i(—ni + h)X 7 
+ |(1 + k - m 3 )I + [k 3 + \(ko + m 3 - 1)] A 3 
+ j-m 3 + |(1 - m 3 - fco)] A' 8 + i(l - ko + m 3 )A' 15 . 



(6.9) 
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In the diagonal part of (|6.9p . there are four independent matrices because equation ()6.9 
represents 4x4 matrix with the structure 



G 



GL{3,C) 
1 



To deal with the matrices from GL(3,C), in the diagonal part of (16. 9p one should separate 
only a 3 x 3 block: 



1)] Af 



Diag = |(1 + k - m 3 )/ (3) + [k 3 + l(ko + m 3 - i;j ^ 3 

+ [-m 3 + |(1 - ma - fed)] A' 8 (3) + ±(1 - fc + m 3 )A; ( 5 3) 

1 
1 
1 



= |(1 + k - m 3 ) 
+ [-m 3 + |(1 - m 3 - fco)] 



+ [A; 3 + i(fco + m 3 -l)] 

+ i(l-fc + m 3 ) 




-1 
1 




-1 



0-10 




Resolving X'jf* in terms of /( 3 ) , a! 3 ^ , A 



'(3). 



A 



/(3) 
15 



we arrive at a 3-term relation: 

Diag = |(1 + 2k - 2m 3 )/ (3) + [A* + ±(fco + 2m 3 - 1)] A^ 3) + §(-4m 3 - 2k + 2)A' £ 

The group law for parameters of SL(3, C) has the form (the notation M = 1 — 2m 3 is used) 

fc ' = k' k + k'k + 2(— n' lo + n'^i), 

(k")i = (k' Q k + k'k + ifc' x k)i + 2(-re Zi + n[l ), 

(k") 2 = (k' k + k'k + ife 1 x fc) 2 + 2(-in' h - in'Jo), 

(k") 3 = (k' Q k + k'k + ik! x k) 3 + 2(-n' l - n^h), 

n o = ( k 'o + k' 3 )n - (k[ - ik' 2 )m + n' M, 

n'{ = (k' Q - feg)ni - (fei + ik' 2 )n + n^M, 

# = l' (ko + fc 3 ) + l[(h + ik 2 ) + M'/ 0) 

l'( = l' (h -ik 2 ) + l[(k -k 3 )+ M'h, 

M" = M'M - 4(Z n - ^m). 

These rules determine multiplication of the matrices 



G 



k + k 3 
k\ + ik 2 
-2l 




k\ — ik 2 
k - k 3 
-2h 




+2n 

-2ni 

M 

1 



If additionally, in equation (16.81) one requires 

tiq = 0, n\ = 0, = 0, Zi = 0, m 3 = 0, uiq = 1, 
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then 



G 



ko + ^3 ki — ik 2 

k± + ik2 +ko — k 3 

10 

1 



with the decomposition rule 

G = fciAi + k 2 X 2 + i(l + ko)I + [k 3 - |(1 - fto)] As + |(1 - fc )A / 8 + |(1 - A ;o )A / 15 - 
One can readily verify that the 2x2 block is given by 

G^{k a ) = hoi® + h\® + k 2 X (2) + k 3 xf\ 

7 On the multiplication law for GL(4, C) in Dirac basis 

In the Gell-Mann basis Aj, an element of GL(4, C) is 

G = a I + (ax - ibt)\i + (o2 - * 2 )A 2 + («A) - 5 3 )A 4 + (A 3 + iBo)A 5 + (A 2 - Bx)X 6 
+ (A x + S 2 )A 7 + (a 3 - ib 3 - ib )X 3 + (-ib 3 + ib - a 3 )A 8 + (-B 1 - A 2 )X g 
+ (Ax - B 2 )X 10 + (iA + S 3 )An + (-A 3 + iB )X 12 + (-ax - ih)X 13 
+ (-a 2 - ib 2 )Xx4 + (ib + a 3 + i& 3 )A' 15 , 

or in variables (k, m, I, n): 

G = i(fco + m )/+ fciAi + /c 2 A 2 + \{(n Q - n 3 ) - (l + Z 3 )]A 4 + ^[-(no - n 3 ) - (*o + Zs)]A 5 
+ |[-(ni + m 2 ) - (Zi - U 2 )]X 6 + ^[(ni + m 2 ) - (Zi - iZ 2 )]A 7 + [k 3 + ~(&o - m )] A 3 
+ [-m 3 + |(m - fco)] A' 8 + |[-(ni - in 2 ) - (Zi + U 2 )]X 9 
+ ^K n i ~ in z) ~ + ^ 2 )]Aio + § [(n + n 3 ) - (Z - Z 3 )]An 
+ 2l[~( n o + "-3) - (Zo - Z 3 )]Ai 2 - miAi 3 - m 2 Ai4 + [m 3 + |(m - k )] X' 15 . 

The problem is to establish the multiplication rule G" = G'G in A-basis: 



As by definition the relationships A m A n = e mn kXk must hold, the multiplication rule is 

The main claim is that the all properties of the GL(4, C) with all its subgroups are determined 
by the bilinear function (|7.1h . the latter is described by structure constants e mn k- It is evident 
that these group constants should be simpler in the Dirac basis Aj than in the basis Aj. Our 
next task is to establish the multiplication law G" = G'G in A-basis: 

-A m A-n — E mn i~A-ki Xfc — E mn f t X m X n . 

Before searching for structural constants E mn ] t , let us introduce a special way to list the Dirac 
basis Aj: 

o>x = 7°7 2 , a 2 = i7°7 5 , a 3 = j^j 2 a 2 = I, axa 2 = ia 3 , a 2 ax = —iax, 




Pi = «7 3 7 1 



= *7 3 



3 = ry 



1 



01 



= I 



(3x(3 2 = i(5 3 ,(i 2 (3x 



03 
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these two set commute with each others ajfik = fyaij, and their multiplications provides us with 
9 remaining basis elements of {A^}: 





= a x 0x = -7 5 , 


Bi 


= a\fh 


= tV, 


Ci 


= ai/? 3 = 7 3 7 5 , 


A 2 


= a 2 f3\ = -ij 2 , 


B 2 


= a 2 02 


= -il 1 ! 2 , 


c 2 


= a 2 3 = -«7 2 7 3 


A 3 


= a 3 0x = 7°, 


B 3 


= a 3 (3 2 


= 7V, 


c 3 


= «3/?3 = 7°7 3 - 



The multiplication rules for basic elements 



ax, a 2 , a 3 


01, 02, 03 




A 


1 


A 2 , A 3 , 


Bx, B 2 , B 3 , 






Cx, C2, C 3 , 






ax 


a 2 


a 3 








01 02 03 




Ax 




A 2 


A3 




ax 


I 


ia 3 


—ia2 




ot-x 




Ax Bx Cx 


ax 


0i 






iA 3 


-iA 2 




a 2 


-ia 3 


I 


iax 








A 2 B 2 C2 


a 2 


-iA 3 




0i 


iAx 




a 3 


ia 2 


—iax 


I 




«3 




A 3 B 3 C 3 


0:3 


iA 


2 




-iAx 


0i 






Bx 


B 2 


B 3 








Cx C2 


c 3 








ax 


a 2 a 3 




ax 


02 


iB 3 


-iB 2 




ax 


03 iC 3 


-iC 2 









Ax 


A 2 A 3 




a 2 


-iB 3 


02 


iBx 




a2 


-iC 3 3 


iCx 







2 


Bx 


B2 B 3 




03 


iB 2 


-iBx 


02 




a 3 


iC 2 -iCx 


03 







\ 


Cx 


C2 c 3 






01 


02 


03 








Ax A 2 


A3 








Bx 


B 2 


B 3 


01 


I 


103 


-i02 


01 




ax Q?2 


a 3 





i 




iCx 


iC 2 


iC 3 


02 


-103 


I 


10X 


02 




-iCx -iC 2 


-iC 3 





2 




O-X 


«2 


0L3 


03 


i02 


-i0x 


I 


03 




iBx iB 2 


iB 3 









-iAx 


-iA 2 


-iA 3 




Cx 


c 2 


c 3 








ax a2 


a 3 








0i 


02 


03 


01 


-iBx 


-iB 2 


-iB 3 




A 


L 


01 iA 3 


-iA 2 




Ax 


ax 


iCx - 


iBx 


02 


iAx 


iA 2 


iA 3 




A 2 


-iA 3 Px 


iAx 




A 2 


a 2 


iC 2 - 


iB 2 


03 


ax 


a 2 


a 3 




A3 


iA 2 -iAx 


0i 




A, 


a 3 


iC 3 - 


%B 3 




Ax 


A 2 


A 3 








Bx B 2 


B3 






Cx 


c 2 


c 3 


Ax 


I 


ia 3 


—ia2 




Ax 


103 -C3 


c 2 


Ax 






102 


B 3 - 


-B 2 


A 2 


-ia 3 


I 


iax 




A- : 




C 3 10 3 


-Cx 


A 2 






B 3 


-i02 


Bx 


A3 


ia 2 


—iax 


I 




-4, 




—C2 Cx 


103 


A, 




B 2 


-Bx 


-i02 




ax 


a 2 


a 3 








01 02 


03 






A 


1 


A 2 A 3 


Bx 


02 


iB 3 


-iB 2 




E 


1 


—iCx ax 


iAx 


Bx 




-103 


c 3 - 


c 2 


B 2 


-iB 3 


02 


iBx 




B 


2 


-iC 2 a 2 


iA 2 


B 2 




-C3 


-i03 Cx 


B3 


iB 2 


-iBx 


02 




B 


3 


-iC 3 a 3 


iA 3 


B 3 




c 


2 


-Cx - 


103 




Bx 


B 2 


B 3 








Cx C 2 


c 3 








ax 


«2 


a 3 


Bx 


I 


ia 3 


—ia2 




Bx 


iPx ~A 3 


A 2 


Cx 




03 


iC 3 


-iC 2 


B 2 


-ia 3 


I 


iax 




B, 




A 3 i0x 


-Ax 


c 2 




iC 3 


03 


iCx 


B3 


ia 2 


—iax 


I 




B3 


-A 2 Ax 


i0i 


c 3 


iC 2 


-id 


03 




01 


02 


03 




Ax A 2 A 3 




B 


1 


B 2 


B 3 


Cx 


iBx 


-iAx 


ax 


Cx 


i(5 2 -B 3 B 2 


Cx 




10X 


A 3 


-A 2 0x 


c 2 


iB 2 


-iA 2 


«2 


c 2 


B 3 i0 2 -B 


1 c 2 




A 3 


-%0x Ax 


C3 


iB 3 


-iA 3 


«3 


c 3 




B 2 Bx i0 2 


c 3 


A 


2 


-A 


x ~i0i 
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Ci 


c 2 


C 3 


Ci 


I 


ia 3 


—ia 2 


c 2 


-ia 3 


I 


ia x 


c 3 


lOL 2 


—ia\ 


I 



(7.3) 



These relations provide us with simple formulas for fifteen coordinates of the element of GL(4, C) 

G = 7 I + ajaj + bjfy + XjAj + YjBj + ZjCj, 

7 = |SpG, a,j = \SpajG, bj = \Sp/3jG, 

X j = \SpA j G, Y j = \SpB j G, Zj = \Sp CjG. 

With the use of relations (|7.3p an explicit form of the group law for (15 + 1) parameters can be 
found: 

( 7 7 + a^Oi + b^ + X'iJU + YlBi + Z'fi^l + a jaj + bjfy + XjAj + + ZjCj) 

= y 7 j + (y aj - + a; 7 )«i + (A- + ^7)& + ( 7 %- + x^Aj + (y^ + y/ 7 )^ 

+ {j'Zj + Zjj)Cj + a' x {a x + a 2 io;3 - «3«a 2 ) + a 2 (—a x ia 3 + «2 + a 3 zai) 

+ a 3 (oiia2 - a2?ai + 03) + a x (b x A x + b 2 B x + &3C1) + a' 2 (b x A 2 + 62-62 + 63C2) 

+ a' 3 (b x A 3 + &2# 3 + 63^3) + a' x (X x (3 x + x 2 iA 3 ~ x 3iA 2 ) 

+ a' 2 (-X x iA 3 + X 2 /3 X + X 3 iA 1 ) + a' 3 (X x iA 2 - X 2 iA 1 + X 3 /3i) 

+ aUn/32 + >2^ 3 - y 3 iB 2 ) + a' 2 (-Y x iB 3 + Y 2 f3 2 + Y 3 iB x ) 

+ a' 3 (yiiS 2 - y 2 iB! + Y 3 (5 2 ) + ai(Zi/% + Z 2 iC 3 - Z 3 iC 2 ) 

+ a! 2 (-Z x iC 3 + Z 2 f3 3 + Z 3 iCi) + 4(ZiiC 2 - Z 2 iC x + Z 3 /%) 

+ b' x (a x A x + a 2 ^2 + 03^3) + b' 2 (aiBi + a 2 -B 2 + 03-63) + b' 3 (a x C x + a 2 C 2 + a 3 C 3 ) 

+ + b 2 ij3 3 - b 3 i(3 2 ) + b' 2 {-b x ip 3 + b 2 + + b' 3 {b x i(5 2 - b 2 i/3 1 + b 3 ) 

+ b' x (X x a x + A 2 a 2 + ^ 3 a 3 ) + &£(-Xi*Ci - X 2 iC 2 - X 3 *C 3 ) 

+ b' 3 (X x iB x + X 2 iB 2 + X 3 iB 3 ) + &i(^'Ci + y 2 iC 2 + y 3 *c 3 ) 

+ 6 2 (yiai + y 2 a 2 + y 3 a 3 ) + 63 (-YiiAi - Y 2 iA 2 - Y 3 iA 3 ) 

+ b' x {-Z x iB x - Z 2 i-B 2 - Z 3 i-B 3 ) + ^(ZiiAx + Z 2 iA 2 + Z 3 ii4g) 

+ b' 3 (Z x a x + Z 2 «2 + Z 3 a 3 ) + X^ai/Si + a 2 iA 3 - a 3 iyl 2 ) 

+ X 2 (-aiM 3 + a 2 /3i + a 3 iili) + X' 3 ( ai iA 2 - a 2 iA x + a 3 /3 x ) 

+ X[(b x a x + b 2 iCi - b 3 iB x ) + X' 2 (b x a 2 + b 2 iC 2 - b 3 iB 2 )+ X' 3 (b x a 3 + 6 2 iC*3 - b 3 iB 3 ) 

+ X[(X X + X 2 ia 3 - X 3 ia 2 )+ X' 2 (-X x ia 3 + X 2 + X 3 ia x )+ X' 3 (X x ia 2 - X 2 ia x + X 3 ) 

+ X' x (Y x if3 3 - Y 2 C 3 + Y 3 C 2 ) + X' 2 (Y X C 3 + Y 2 i(3 3 -Y 3 C x ) + X' 3 (-Y X C 2 + Y 2 C x +Y 3 i(3 3 ) 

+ X' x (-Z x ip 2 + Z 2 B 3 - Z 3 B 2 ) + X' 2 {-Z X B 3 - Z 2 i(3 2 + Z 3 B X ) 

+ X' 3 (Z X B 2 - Z 2 B X - Z 3 ip 2 ) + Y((a x p 2 + a 2 iB 3 - a 3 iB 2 ) 

+ y 2 (— aiiB 3 + a 2 [3 2 + a 3 iB x ) + y 3 '(aii5 2 - a 2 iB x + a 3 (5 2 ) 

+ Y((-b x iC x + 6 2 "i + 63*^1) + Y 2 (-b 1 iC 2 + b 2 a 2 + b 3 iA 2 ) 

+ y 3 '(-MC3 + & 2 a 3 + 6 3 *40 + Y{{-Xiifa + X 2 C 3 - X 3 C 2 ) 

+ Y 2 '(-X X C 3 - X 2 ifa + X 3 C X ) + Y 3 '(X X C 2 - X 2 C X - X 3 ifo) 

+ Y((Y X + Y 2 ia 3 - Y 3 ia 2 ) + Yi{-Y x ia 3 + Y 2 + Y 3 ia x ) + Y 3 \Y x ia 2 - Y 2 ia x + Y 3 ) 

+ Y{{Z x i(5 x - Z 2 A 3 + Z 3 A 2 ) + Yi{Z x A 3 + Z 2 ij3 x - Z 3 A X ) 

+ Y 3 \-Z X A 2 + Z 2 A X + Z 3 i(3 x ) + Z' x {a x l3 3 + a 2 iC 3 - a 3 iC 2 ) 

+ Z 2 (-aiiC 3 + a 2 j3 3 + a 3 iC x ) + Z 3 (aiiC 2 - a 2 iC x + a 3 (5 3 ) 
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+ Z' x (b x iB x - b 2 iA 1 + 6 3 ai) + Z 2 (b x iB 2 - b 2 iA 2 + b 3 a 2 ) + Z' 3 (b x iB 3 - b 2 iA 3 + b 3 a 3 ) 
+ Z' x {X x if3 2 - X 2 B 3 + X 3 B 2 ) + Z^Xiflg + *2*#2 - X 3 B X ) 
+ Z' 3 {-X X B 2 + X 2 B X + X 3 i/? 2 ) + Zii-Yrffc + y 2 ^ 3 - y 3 ^4 2 ) 
+ Z' 2 (-Y X A 3 - Y 2 ip x + y 3 ^i) + Z' 3 (Y X A 2 - Y 2 A X - Y 3 ifc) 

+ Z[(Z X + Z 2 ia 3 - Z 3 ia 2 ) + Z 2 (-Z x ia 3 + Z 2 + Z 3 ia x ) + Z 3 (Ziia 2 - Z 2 iai + Z 3 ). 
From these relations we arrive at the following composition rules: 

7" = il + Kai + a' 2 a 2 + a' 3 a 3 ) + (6^6 a + b' 2 b 2 + 6' 3 6 3 ) 

+ (X' X X X + X' 2 X 2 + X' 3 X 3 ) + {Y{Y l + Y 2 % + Y3V3) + {ZiZt + Z' 2 Z 2 + Z' 3 Z 3 ), 
a'[ = {ia x + a' l7 ) + {b' x X x + ^1 + W + + Y[b 2 + Z(6 3 ) 

+ i(a' 2 a 3 - a' 3 a 2 ) + i(X^X 3 - X' 3 X 2 ) + *(y 2 'y 3 - Y 3 "Y 2 ) + i(Z' 2 Z 3 - Z' 3 Z 2 ), 
a 2 = (l'a 2 + a' 27 ) + {b\X 2 + b' 2 Y 2 + b' 3 Z 2 ) + (X' 2 b x + y 2 6 2 + Z' 2 b 3 ) 

+ i(a' 3 a x - a' ia3 ) + i(X 3 Xx - X;X 3 ) + i{Y 3 'Y x - Y{Y 3 ) + i{Z' 3 Z x - Z[Z 3 ), 
4 = ( 7 'a 3 + a' 37 ) + + 6' 2 y 3 + 6' 3 Z 3 ) + (X' 3 b x + Y 3 % + Z 3 6 3 ) 

+ i(a' ia2 - a' 2ai ) + i(X[X 2 - X' 2 X X ) + i(Y{Y 2 - Y 2 'Y X ) + i{Z[Z 2 - Z' 2 Z X ), 
b'{ = 761 + &17 + i(6' 2 fe 3 - 6' 3 fe 2 ) + (aiX x + a' 2 X 2 + a' 3 X 3 ) + (X[ai + X 2 a 2 + X 3 a 3 ) 

+ i(Y(z x + y 2 % + y 3 %) - i(z;y + z' 2 y 2 + z 3 y 3 ), 

6 2 = l'b 2 + 6 27 + i(b' 3 h - b\b 3 ) + (a^y + a' 2 Y 2 + a 3 y 3 ) + (Y[a x + y 2 'a 2 + y 3 'a 3 ) 

+ i(Z' x X x + Z 2 X 2 + Z' 3 X 3 ) - i(X[Z x + X 2 Z 2 + X' 3 Z 3 ), 
b 3 = j'b 3 + 763 + i(b' x b 2 - &2&1) + (a'l^i + a 2 Z 2 + a' 3 Z 3 ) + (Z(ai + Z 2 a 2 + Z 3 a 3 ) 

+ i(XiYi + X' 2 Y 2 + X' 3 Y 3 ) - i(Y(X x + Y 2 'X 2 + ypf 3 ), 
X'( = ( 7 'X! + 7XO + (a'ifti + ai6i) + z(y/6 3 - y^) + - 6 2 Z[) 

+ i(a' 2 X 3 - a' 3 X 2 ) - i(a 2 X' 3 - a 3 X' 2 ) + (Z 2 Y 3 ' - Z 3 Y 2 ') + (Z' 2 Y 3 - Z' 3 Y 2 ), 
X' 2 ' = (iX 2 + 7 X' 2 ) + + a 2 6;) + i(y 2 '6 3 - Y 2 b' 3 ) + z(6' 2 Z 2 - 6 2 Z 2 ) 

+ i(a 3 X x - a;X 2 ) - i(a 3 X; - ai X 3 ) + (Z 3 Y? - Z{4) + [Z' 3 Y X - Z[Y 3 ), 
X'i = (j'X 3 + 7 X 3 ) + (a' 3 h + a^) + i(Y 3 % - Y 3 b' 3 ) + i(b' 2 Z 3 - b 2 Z' 3 ) 

+ i(a! x X 2 - d 2 X x ) - i( ai X' 2 - a 2 X[) + (Z X Y 2 ' - Z 2 Y[) + (Z[Y 2 - Z' 2 Y X \ 

y/' = (yy + 7 y/) + (a[b 2 + a x b' 2 ) + i(z' x b x - M) + i{b' 3 x x - b 3 x[) 

+ i(a' 2 Y 3 - a' 3 Y 2 ) - i(a 2 Y 3 ' - a 3 Y^) + (X 2 Z' 3 - X 3 Z' 2 ) + (X' 2 Z 3 - X' 3 Z 2 ), 

y 2 = (1% + 7I2) + (a' 2 &2 + « 2 fo 2 ) + *(^ 2 &i - + i(b' 3 X 2 - b 3 X' 2 ) 

+ i{a! 3 Y x - a' x Y 3 ) - i{a 3 Y{ - a x Y^) + (X 3 Z' X - X X Z' 3 ) + (X' 3 Z X - X' X Z 3 ), 
y 3 " = (7% + 7^3) + (4&2 + a 3 6 2 ) + i(Z' 3 b x - Z 3 b' x ) + i(6' 3 X 3 - b 3 X' 3 ) 

+ i(a' i y 2 - a' 2 Y x ) - i{a x Y 2 ' - a 2 Y() + [X X Z' 2 - X 2 Z[) + {X[Z 2 - X' 2 Z X ), 
Z'{ = (7% + 1 Z' X ) + (a\b 3 + a x b' 3 ) + i(Y x b' x - Y{b x ) + i(X[b 2 - X x b' 2 ) 

+ i(a' 2 Z 3 - a' 3 Z 2 ) - i(a 2 Z' 3 - a 3 Z' 2 ) + (Y 2 X' 3 - Y 3 X' 2 ) + (Y 2 'X 3 - Y 3 'X 2 ), 
Z'l = {iZ 2 + 1 Z' 2 ) + (a' 2 b 3 + a 2 b' 3 ) + i(Y 2 b' x - Y 2 % x ) + i(X' 2 b 2 - X 2 b' 2 ) 

+ i(a' 3 Z x - a' x Z 3 ) - i(a 3 Z' x - a x Z' 3 ) + (Y 3 X[ - Y X X' 3 ) + (Y 3 'X X - Y(X 3 ), 
Z'i = (7% + 1 Z' 3 ) + (a' 3 b 3 + a 3 b' 3 ) + i(Y 3 b' x - Y 3 'b x ) + i{X' 3 b 2 - X 3 b' 2 ) 

+ i(a' x Z 2 - a' 2 Z x ) - i(a x Z' 2 - a 2 Z[) + (Y x X f 2 - Y 2 X[) + (Y{X 2 - Y^X X ). (7.4) 

With the help of the index notation 

X = CW, Y = C^ 2 \ Z = C$\ 
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it is easy to see a cyclic symmetry in the above relationships: 

7 --f-f + a k a k + b k b k + C k C k + L k L k +L k L k , 

4 = ia k + 7 a' fc + (&;C7« + hcl iy ) + (b' 2 C ( k 2) + b 2 cf) + (b' 3 C { k 3) ) + 6 3 cf' ) 

+ ie kln a' ian + ie kln ci lY C^ + ie Hn c\ 2) ' + ie kln cf Y C®, 
K = l'bk + lb' k + ie Hn b\b n + Kcf^ + ai c[ kY ) + (a 2 C? } + a 2 C? } ') 

+ (a 3 cf } + a 3 C? } ') + ieu n C®'C&\ 
Cf = T'^ 1} + 7^ 1} ' + (a' fc 6i + a fc &;) + ie (1)In (cf 6„ - cf &'J 

+ ie^C^ - a,C«') + e kln {c[ 2) ' + c\ 2) C^'), 
Cf = 7'cf + 7Cf' + (a'fcba + a k b' 2 ) + ie (2) , n (cf '&„ - cf b' n ) 

CP ~ mCW) + e kln (cfc^ + C®CW), 
Cf = 7 'cf } + 7 cf' + (46 3 + afcb's) + ie(3)in(C { k lY b n - C®b' n ) 

+ iekMC^ - ai C^') + e kln (cl lY C^ + C^C^'). (7.5) 

It is readily seen that these group multiplication laws (|7.4p . (|7.5[) permit 15 two-parametric 
subgroups: 

(7,0) G {(7,01), (7,a 2 ), (7,03), (7, 61), (7, 6 2 ), 7> &3), (7, -^1), (7> X 2 ), (7, X 3 ), 
(7, Yi), (7, Y 2 ), (7, 13), (7, Zi), (7, Z 2 ), (7, ^3)} 

with the same composition law: 

l" = 7 7 + a ' a > a " = 7 a + 7 a 'i 

which in variables 7 = W cos 4>, a = iW sin takes the form 

VP" = W'W, a" = a + a. 

The variable W is determined by det G(W, a) = W 4 , the choice W = 1 guarantees det G = +1. 
All 15 basis elements A( p ) G {a k , (3 k , A k , B k ,C k } possess the same properties: 

A Co)= A (p)> A (p) =/ - 

Therefore, one can construct 15 different elementary unitary (at real valued parameters) matrices 
by one the same recipe: 

U ( P ) = e**W A W = cos^ (/3) +isin^ (p) A (p) , 

K) = U (p) = e ~ lK)hp) = COS< ^) " ^m^ (p) A (p) . 

The whole set of unitary matrices SU{4) may be constructed on the basis of a simple factorized 
formula: 

U = e ^(i) A (i) . . . e ^d5) A (i5)_ 



The order of the factors is important. Every such order leads us to a definite parametrization 
for the group SU(4) - all them seem to be equivalent. 
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In the end of the section let us write down the explicit form of these 15 elementary unitary 
transformations: 



Oil 



An 



Bo 



Co 



a 2 


COS ( 





~(J 2 

i + i sin <fi(j2 




«2 







COS ( 



I 

i 


i sin <J)(j2 



a 3 



a 2 
o 2 



COS ( 

sin ( 



U 



cos cp I sin (pG2 
i sin (f)(j2 cos (p 

—io\ 
io\ 

cos (j) sm (pas 
— sin (pas cos (p 

ia 2 
ia 2 

cos (f> — sm <poi 
sin 0(72 cos <p 



-cr 3 

- i sin (pas 




/9 



COS i 



a 2 
a 2 

+ i sin <pa2 




02 



— sm 
cos q 

- 
icr 3 



ias 








-4, 



/ 
/ 



COS <p + I Sin 0(72 

cos <p sin 0(7i 
sin <pa\ cos 

cos 



cos 9 
zsini 



? sm ( 
cos a 



+ 1 sm i 


Si 



cos <p 


0"! 



jsini 



0"! 





-cr 3 


cos c 



B, 



COS i 



-(71 



COS ( 



- 1 sm <pa\ 




-(71 

- i sin 0o"i 




cos (p + i sin 0(73 




COS i 
COS i 

Cs = 

COS ( 
COS ( 



-(71 
(71 



? sin (pas 


i sin 0<7i 


-cr 3 



i sin <pa\ 


% sin <^><7 3 



U 



B 



cos i sm 0o"i 
i sin <pa\ cos 



CT3 




-cr 3 






-cr 3 

COS 

— i sin 0(73 



-J Sill 0(73 
COS (p 



(7.6) 



Certainly, these relations provide us with 15 elementary solutions of the unitarity equations (|3.3p . 
For instance, the generator 0:2 gives rise to the above 1-parametric Abelian subgroup G°(a); 
whereas the above 4-parametric subgroup Go x SU(2) (|5.15p is generated by (02; 0i, B2, C2). 

The question is how one could describe all combinations of the above 15 simple sub-solutions 
by a single unifying formula - the latter should evidently exist. 



8 On factorization SU (4) and the group fine-structure 

On the basis of 9 matrices (I7.2p one can construct six 3-dimensional sub-sets: 

K = {A 1 = a 1 (3 1 , B 2 = a 2 /?2, C 3 = a 3 (3 3 }, 

L = {d = aifls, A 2 = a 2 fli, B 3 = a 3 (3 2 }, 

M = {B\ = anfo, C 2 = a 2 p 3 , A 3 = a 3 /9i}, 

K' = {-C[ = -otifo, -B' 2 = -a 2 /?2, -C's = -a 3 /3 3 }, 

L' = {-B[ = -ai/3 2 , -A' 2 = -a 2 Pi, -B' 3 = -a 3 /? 2 }, 

M' = {-A', = -ai/3i, -C' 2 = -a 2 /%, -B' 3 = -a 3 /3 2 }, 
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(one may recall the rule to calculate the determinant of a 3 x 3 matrix) with the same commu- 
tation relations: 

rir 2 = -r 3 , r 2 r 1 = -r 3 , r^-r^^o, + r 2 ri = -2r 3 , (8.1) 

and analogous by cyclic symmetry. The whole set of the above 9 matrices coincides with 

a, (3, K, L, M, (8.2) 

or 

a, 0, K', L', M'. 

It suffices to consider one variant, let it be (|8.2p . It seem reasonable to suppose that arbitrary 
element from GL(4, C) can be factorized as follows 

S = e im e*$e ikK e <lI 'e imM . (8.3) 

When all parameters are real-valued, the formula provides us with the rule to construct elements 
from SU(4) group0. The order of factors might be different. Let us specify the group law for 
these 5 subsets. First are the two groups: 

e ma = cos a + i sina(niQi + n2«2 + n^a^), e lbl3 = cos b + i sin b(mf3\ + ^2/^2 + nzfls)- 
They are isomorphic, so one can consider only the first one: 

e ma = cos a + isina(niai + 71202 + ^303) = %o — ix\a\ — 1x20-2 — ^as- 
Multiplying two matrices we arrive at 

Xq = XqXo — x\x\ — x' 2 X2 — X3X3, x'l = XqXi + X^Xq + (x' 2 Xs — X2X3), 

A = x 'o x 2 + x' 2 xq + (x' 3 xi - x[x 3 ), x'l = x' x 3 + x' 3 x + (x[x2 - x' 2 xi). (8.4) 
Parameters (xo,Xi) should obey 

xl + x\ + x\ + x\ = 1 <=^ det e iS " = +1. 
The inverse matrix looks 

{x Q ,x)~ l = (x ,-x). 

With real (xo,Xi) we have a group isomorphic to SU(2), spinor covering for SO(3,R): 

x .. d + c + c' x c 

c=— , c = . 

Xo 1 — c'c 

At complex (xo,Xi) we have a group isomorphic to GL(2,C), spinor covering for SO(3,C) or 
Lorentz group. 

Now let us turn to finite transformations from remaining subsets. It is readily verified that 
these 1-parametric finite elements 

e «yiri _ COS y 1 _|_ isinyiTi, e i v^i = COS y 2 -f isaxy^^i e mT3 = cosy^ + ismy 3 T 3 , 



Just such a structure was described in |40] 
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commute with each other: 



e tyiVl e iy2V2 = (cosyi + isinyiri)(cosy2 + ismy 2 Y 2 ) 



= cos yi cos y 2 + i cos yi sin 7/2^ + i cos 2/2 sin yiTi + sin y\ sin y2r3i 
e iy2T2 e iyiTl = (cos 2/2 + isin^I^Xcosyi + i sin yiTi) = 

= cos ?/2 cos yi + i cos ?/2 sin yiTi + z cos j/i sin y 2 T 2 + sin 1/2 sin y\T 3 , 

that is e iyiVl e iy2T2 = e^ 2r2 e l2/iri , and so on. Evidently, this property correlates with the com- 
mutative relations (|8.ip . Thus, each of tree subgroups can be constructed as multiplying of 
elementary 1-parametric commuting transformations. Their explicit forms are: 
subgroup K 



K = {A x = ai/h, B 2 = a 2 f3 2 , C 3 = a 3 (5 3 ], 



ik\K\ 



I 
-/ 



Bo 



cos k\ + i sin feiAi, 
3 ifc 3 ^3 =cosA:3 + isinfc3C3; 



-0-3 






-0-3 



0-3 






-0-3 



ifc2 K\ 



cos k2 + i sin fc 2 -B 2) 



subgroup L 



{d = ai/3 3 , A 2 = a 2 f3i, B 3 = a 3 /3 2 }, 





-0-3 



-0-3 





A, 



cos/i + isin/iCi, 



2<7 2 

-icr 2 

= cosZ 2 + isini 2 A 2 , 



-0-1 

CTl 

e ihL3 _ 



cos Z 3 + i sin /3-B3; 



subgroup M 



M 

Bi 



{B\ = ai(3 2 , C 2 = a 2 /?3, A 3 = a 3 (3i}, 

(71 
CTl 



im\M\ 



cos m\ + i sin mi £>i , 
cos 7713 + i sin 7713^3 . 



-CTl 

-0-1 

im,2 Mi 



A, 



/ 
/ 



cos m 2 + « sin m,2C2 



One additional note should be made. In the recent paper by A. Gsponer [35] on the quater- 
nion approach to the problem of building the finite transformations from SU (3) and SU{4) an 
important point was to divide 15 basis 4x4 matrices into three sets: 

set A of antisymmetrical matrices, 
set S of symmetrical matrices, 
set D of diagonal traceless ones. 

It is easily seen that 

set A = {a { © Pi} ; 

set S = {A 2 , A 3 , B u B 3 ,d,C 2 } = {L ® M}; 
set D = {Ai,B 2 ,C 3 = K}. 
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Turning again to relationship (|8.3p . let us rewrite it as follows 
S 



which exactly corresponds to the structure used in [35] in connection with the Lanczos decom- 
position theorem [43j . 

Several last comments should be made. On the basis of 15 matrices 



Oil, 

A 2 
A 3 



Oi 2 , 
ai fa, 
a 2 Pi, 
a 3 Pi, 



«3, 

Bi 
B 2 
B 3 



Pi, 
aiP 2 , 
a 2 p 2 , 

OL?,P 2 , 



Pi, Pz, 
Ci = aiPz, 
C 2 = a 2 p 3 , 
C 3 = a 3 p 3 



one can easily see the following 20 ways to separate SU{2) subgroups (certainly, those arise at 
real-valued parameters; complex-valued parameters give rise to linear subgroups GL(2,C)): 



(ai, a 2 , a 3 ), 
(a 1 ,A 2 ,A 3 ), 
(a 1 ,B 2 ,B 3 ), 
(ai,C 2 ,C 3 ), 
(Pi,B\,Ci), 
(Ai,p 2 ,Ci), 
(A^B^ps), 



(Pl,P2,P 3 ), 

(Ai,a 2 ,A 3 ), 

(B 1 ,a 2 ,B 3 ), 

(Ci,a 2 ,C 3 ), 

(Pi, B 2 , C 2 ), 

(A 2 ,p 2 ,C 2 ), 

(A 2 ,B 2 ,p 3 ), 



(Ai,A 2 ,a 3 ), 

(B 1 ,B 2 ,a 3 ), 

(Ci,C 2 ,a 3 ), 

(Pi,B 3 , C 3 ), 

(A 3 ,p 2 ,C 3 ), 

(A 3 ,B 3 ,p 3 ). 



(8.5) 



Certainly, they provide us with twenty different 3-parametric solutions of the unitarity equa- 
tions (|3.3p . Such 3-subgroups might be used as bigger elementary blocks in constructing of a 
general transformation [25, 28]. 

For instance, for the variant from (|8.5p : (a\, A 2 , A 3 ) =^ ( a i, X 2 , X 3 ) the general multiplica- 
tion law (|7.4|) gives 



/ = 7 ' 7 + a[ ai + X' 2 X 2 + X' 3 X 3 , a!{ = j'ai + o' l7 + i(X' 2 X 3 - X' 3 X 2 ), 

a' 2 ' = 0, 4 = 0, b'( = 0, b% = 0, b' 3 ' = 0, X[' = 0, 

X'J, = iX 2 + 7 A 2 + +i{-a! x X 2 ) + aiX' 3 ), X' 3 ' = 7 ' ' X 3 + 7 X' 3 + i(a! x X 2 - ai X' 2 ), 

y(' = o, y 2 " = o, y 3 " = o, z? = o, z» = o, z 3 ' = o, 

that is 

7" = 77 + a^ai + X' 2 X 2 + X 3 X 3 , a" = 701 + a' l7 + i(X' 2 X 3 - X 3 X 2 ), 
X'{ = 7 'A 2 + 7 X' 2 + i(X' 3 ai - aiX 2 ), X' 3 ' = 7 'X 3 + 7X3 + ifaXz - aiX' 2 ), 

which coincides with equation (|8.4p . The same can be done for any other representative 
from (1831). 



9 On pseudo-unitary group SU(2, 2) 

As said, the Dirac basis was used previously [9j |4T)1 l4~8j [50] in studying the exponentials for 
SU(2, 2) matrices. Let us show how the above formalism can apply to this pseudo-unitary 
group SU(2, 2). Transformations from SU(2, 2) should leave invariant the following form 



( * * * * \ 
\z x ,z 2 ,z 3 ,z 4 ) 



+1 

+1 


















Zl 

z 3 

Zl 



z r\z = z T]z , 
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which leads to 



Uz, 



where U + r] = r/U , 







+1 











1 

; 







+1 








v = 








-1 


















-1 



Any generator A' of those transformations must obey relation 

U k = e iaA ^ (A' k )+V = vKi €{!,..., 15} 

or allowing for identity T] = — 7 s 

(A' fc ) + 7 5 = 7 5 A' fc , fce {!,..., 15}. 

All generators A' k of the group SU(2, 2) can be readily constructed on the basis of the known 
generators of SU(4) (see (|6.1j) ) 



a; 

A', 



Ai = r 

iA A = - 



a; = zAi 



Z7 



7 1 , A' 5 = iA 5 = i 7 V, 



Ar 



7 



iAy = ?7 5 7 2 , 



l 5 

A' 



A' 3 = iA 3 
A' = 



A' 10 = A10 = 2a 01 , A' n = A n = 2a 02 , 
l 13 - ,t 13 = 2ia 12 , A' u = A i4 = 2ia 23 , 

Basis elements may be listed as follows: 



A', 



Ai 



iA 8 = 
: A n 
A' = Ai 



-7 



-7 5 7°, 



A' 



zA 6 
iA 9 



-7 , 
i 7 5 7 3 , 



A' 10 = A12 = 2a 



OH 



L 12 

A' 



15 



A 



15 



2ia 



31 



(9.1) 



Q 1 = OL\ 

K) 2 = - 

a' 2 a' 3 = - 



7V, 



a 2 = 1012 



I \2 



K) 



-7°7 5 , 

' \2 _ 



a 3 = za 3 



Z7 5 7 2 



K) 



3"1 



za 



2' 



a, a 



1"2 



za 



3- 



# = A = irV, /3 2 = i/3 2 = -7 3 , /3 3 = ^3 = -7 1 , 



1/ \2 



(ft) 



(ft) 



/ \2 



(ft) 



These two sets commute with each other: ctjfik = fto^; an d their multiplications give remaining 
9 elements: 



iA 2 



1 o> 



"2ft = 7 



4 = iA 3 = ol 3 x = 17 







51 



-B 2 = 
-fl 3 



"lft = 

-- olA 



+Z7V, 

= -7°y, 



C[ = id = a[f3' 3 



i 7 3 7 5 , 



C 2 



-Co 



«2P3 



■C3 = ^ 



+i 7 2 7 3 , 
-7°7 3 - 



Making in relations (|7.6p a formal change in accordance with 
e mA = cos a + z'sinaA, A' = zA, 

cos a + i sin a A = e iaA ==>• e iaA = cos z'a + i sin z'a A = cosh a — sinh aA, 



we arrive at explicit form of elementary pseudo-unitary SU(2, 2)-transformations: 



/ 

ct 1 



(To 





tt 2 = ZQ; 2 





a 2 

-1 

1 



cos (j) + i sin <pa2 

cos 1 





- i sin (/)0"2 



coshx — isinhx 
i sinh x cosh \ 
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a 3 = ia 3 
&=Pl = 
02 = ifa ■ 
0' 3 = i/% ■ 

A[ = A, : 

A'o = iA 2 



A' 3 = iA 3 



B' 



B> 



-B 2 
-B 3 



C[ = id = 

c' 2 = — c% 



-c 3 



ia 2 
ia 2 

(72 

a 2 

(7 3 
-a 3 

(71 
-(71 

/ 
-/ 

-(7 2 
(7 2 

i 

1 

io\ 
io\ 

(73 
(7 3 

(7i 
-(7i 



-ia 3 

(71 
(7i 

"(73 
+(7 3 



U: 



B 



ia 3 




cosh x — sinh X°2 
sinh XC2 cosh x 

cos + i sin 0(72 

cos ^ + i sin 002 

cosh x i sinh XC3 
z sinh \a 3 cosh x 

cosh x i sinh x^i 
z sinh x&i cosh x 

cos (f) + i sin </> 

cos (f> — i sin 

coshx — isinhx(72 
z sinh x(72 cosh % 

cosh x ~ sinh x 
- sinh x cosh x 

cosh x — sinh x(7i 
— sinh X&1 cosh % 

cos 4> + i sin ^(73 

cos + z sin </>o"3 

cos ^ + i sin 0o"i 

cos (f> — i sin <pa± 

cosh x sinh \a 3 
sinh XC3 cosh x 

cos <f>i sin </><7i 

cos 4> + i sin </><7i 



u§ 



cos — i sin ^(73 

cos 4> + i sin ^(73 



(9.2) 



We can easily obtain unitarity equations for SU(2, 2) group, simple solutions to which are 
given by (|9.2p . Indeed, taking into account the formulas (see (|3.ip 



G + r] -- 
r/G- 1 



kl + k*a 



-R - I* a 



11 







n a m — ma 



-I 
/ 



k' + k'a 
— 1' — I' a m\ 



-I 
+1 



n'a 



ma 



-k 







k*a 



Uq + n*a m 



-l* -I* a 
m*a 







-k' — k'a —n' + n'a 
—1' — I' a m' — m'a 



from G + rj = r/G 1 we produce 



h* — k' 



I* 
'0 



0- 

n' , 



0- 



V 



n 



n 



— '0' 



In. 



m 

* 

n - 



m 
-I'. 



These relations differ from analogous ones (|3.2p for SU (4) group only in all signs of the second 
line. Therefore, the unitarity conditions for the group SU(2,2) are (compare with (|3.3p ) 



/cq = +ko(mm) + mo(ln) + lo(nm) — no(lm) + il(m x n), 
m,Q = +mo(tt) + ko(nl) + no(lk) — lo(nk) — in(k x I), 
k* = —k(mm) — m{ln) — l{nm) + n(lm) + 21 x (n x m) 
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+ imo(n x I) + il (n x m) + ino(l x m), 
m* = -m(kk) - k(nl) - n(lk) + l(nk) + 2n x (I x fe) 

— zfco(i x n) — mo(i x k) — Uq(ti x fc), 
— = +A:o(nm) — mo(kn) + Zo( nn ) + no(km) + zfc(n x m), 
— ng = +mo(Ifc) — ko(ml) + tiq(11) + lo(mk) — im(l x fc), 
— Z* = —k{nm) + m(kn) — l(nn) — n(km) + 2k x (m x n) 

+ ik()(m x n) + imo(k x n) + irio(m x k), 
—n* = —m{kl) + k(ml) — n(ll) — l(mk) + 2m x (k x I) — imo(k x I) 

— iko(m x I) — ilo(k x m). 

Several words on factorization SU(2, 2) = SU(1, 1) x [K x L x M] x SU(1, 1) and a further 
group fine-structure for SU(2, 2). On the basis of 9 matrices 

^ = A 1 = ai/3j = - 7 5 , Bi=iBi = ai^ = i 7 5 7 1 > C( = = a'^ = i 7 V , 
A' 2 = z,4 2 = a' 2 # = 7 2 , B 2 = -B 2 = a' 2 f5' 2 = +i 7 V, C 2 = -C 2 = a' 2 {3' 3 = ry 2 7 3 , 
A' 3 = iA 3 = a' 3 # = i 7 °, # 3 = -fl 3 = a 3 /3 2 = - 7 ° 7 1 , C' 3 = -C 3 = a' 3 P' 3 = - 7 ° 7 3 

one can construct three 3-dimensional subsets (omitting three others): 

K' = {A[ =Ai = a'^i, B' 2 = -B 2 = a' 2 f3' 2 , C 3 = -C 3 = a' 3 f3' 3 }, 
L' = {C[ = id = oi^, A' 2 = iA 2 = a' 2 p[, B 3 = -B 3 = a' 3 (5' 2 }, 
M' = {B[ = iB 1 = C' 2 = -C 2 = a' 2 p' 3 , A' 3 = iA 3 = a'M, 

with the same commutation relations: 

rir 2 = +r 3 , r 2 r 1 = +r 3 , rir 2 -r 2 ri = o, + r 2 ri = +2r 3 , 

and analogous ones by cyclic symmetry. Arbitrary element from SU(2, 2) can be factorized as 
follows 

C _ Jaa' Ab/3'ikK'ilL'JmM' 
O — c o c o c- , 

all parameters are real-valued. Let us specify the group law for these 5 sub-sets. Two groups 
e iaa ' , e lb ^' are isomorphic so one can consider only the first one: 

e iact' _ j + j( aiQ ,/ _|_ a2Q ' 2 + a 3 a' 3 ) - 7g(a\- a 2 - a 3 ) 

- ^(a 2 - a| — a|)(aia' 1 + a 2 a 2 + a 3 a 3 ) + ^ (a 2 - a 2 - a§) 2 H . 

In the variables 

2 2 2 i 

dj = anj, n 1 — n 2 — n 3 = 1 
we have 

e MQ = / + i(nia'i + n 2 c/ 2 + n 3 a 3 )a — |,a 2 — ^(nia^ + n 2 a' 2 + n 3 c/ 3 )a 3 
+ h. a + II ( n i a 'i + n 2a' 2 + in 3 a' 3 )a 5 - ^a 6 H , 

that is 

e MQ = cos a + i sin a(nia' 1 + n 2 ol 2 + n 3 c/ 3 ) = xq — ixia^ — ix 2 a 2 — ix 3 a' 3 . 
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Multiplying two matrices we arrive at 

Xq = XqXq - x' x X\ + x' 2 X 2 + X3X3, 

x' x 2 + x' 2 x + (x' 3 x 1 - ^1^3), 



• r 2 



x'( = x' x\ + x[x - (x' 2 x 3 - x 2 x 3 ), 

x 3 = X X 3 + X> 3 X + ( x l x 2 ~ x 2 x l)- 



The inverse matrix looks 

(x ,x)- 1 = (x ,-x). 
Parameters (xQ,Xi) should obey the following condition 
% = l. 



2 1 2 

Xq -f" X^ 



X n 



In the variables o, rij it will look 



cos 2 a + sin 2 a(n 2 



2 2^ 

n — nV 



L 2 ~ '"3J — 1' 
For three particular cases we will have: 



1. 



n = (1,0,0), cos 2 a + sin 2 a = 1, 
n = (0, i, 0), cos 2 a — sin 2 a = 1, 



e wai = cos a + i sin aa^ , a G i?; 



a = i&, 



cos ifr = cosh b, 



s'mib = i sinh 6, b G i?, 



cosh 6 — i sinh 6a' 



2- 



n = (0,0, z), 
a = ib, cosib = coshb 



9 O 

cos a — sin a = 1, 



sin ib = i sinh 6, b £ R, 



e iaa 3 



cos a + i sin 071303 = cosh 6 — i sinh 6a 3 



Now let us turn to finite transformations from remaining sub-sets K' , L' , iW"'. Each of tree 
subgroups can be constructed as multiplying of elementary 1-parametric commuting transfor- 
mations. Their explicit forms are: 

subgroup K' 

K' = {A\ =A X = a[(3[, B 2 = -B 2 = a' 2 (3' 2 , C 3 = -C 3 = a' 3 (3' 3 }, 



1 
—i 



B' 



(73 
CT 3 



subgroup L' 

L' = {C[ = iCi = a[0 3 , A 2 = iA 2 = a 2 x , B' 3 



-a 3 
a 3 



-B 3 = a' 3 P' 2 }, 





-ia 3 



-ia 3 





cr 2 



-cr 2 




B' 



&1 







-CTl 



subgroup M' 



M' = {B[ = iB x = ai fo, C 2 = -C 2 = a 2 p 3 , A' 3 



ia\ 
io\ 



C", 



C7i 
(71 



A' 



iA 3 = a 3 Pi}, 

i 

1 



On the basis of 15 matrices 



"2> 



«3> 



3) 
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4 


/ al 
= "lPl, 




' al 
= OlP2> 




= 4$, 


4 


' 0/ 
= 02Pl> 




' al 

= a 2 P 2 , 


c 2 


1 ai 


4 


= «3Pl) 


4 


1 al 

= a 3P2, 


c 3 


= 4/?3 



one can easily see the following 20 ways to separate SU(l, 1) subgroups: 



(04,04,03), 






(44,4), 




(4,4, 4> 




(^1,4,^3), 


(B[,B 2 ,a 3 


(o:i, C 2 ,C 3 ), 


(C^O^Cg), 


(C[,C' 2 ,a' 3 ) 


((3[, B[,C[), 


(Ai,-^,^), 


(f3'i,B' 3 ,C 3 ) 


(4, P2, Ci)> 


(4, P2, ^2), 


(4> C3) 




(^ 2 , J B 2 ,/3 3 ), 


(^3,-03,03) 



such 3-subgroups might be used as bigger elementary blocks in constructing a general transfor- 
mation. 



10 On pseudo-unitary group SU(3, 1) 



Let us show how the above formalism can apply to the pseudo-unitary group SU(3, 1). Trans- 
formations from SU(3, 1) should leave invariant the following form 



/ * * * * \ 

\ z l , z 2 , z 3 , Z A ) 



+1 













Z\ 





+1 










z 2 








+1 







z 3 











-1 







z rjz = z rjz , 



which leads to 



Uz, 



where U + rj = rjU , 



V 



+ 10 
0+10 
0+1 




Any generator A' of those transformations must obey the relation 

U k = e iaA '*, (A'J+r? = 7]A' k . k G {1, . . . , 15}. 
The matrix rj is a linear combination 

r) = \{2ia 12 - 2a 03 - 7 5 + I) = ±(i 7 V - 7 V + ^SV + T). 

All generators of the group SU (3, 1) can readily be constructed on the basis of the known 
generators A fc of 577(4) (see (|PD ) 



Ai, A 2 , A3, A4, A5, A6, 
iAn, iAi2, «Ai3, iAi4, A15; 



Ay, 



As 



iXg, 



iAio, 



10.1) 



generator A9, . . . , A14 are multiplied by imaginary unit i. Instead of As, A15 one can introduce 
other generator A' 8 , A' 15 see (|6.4p (diagonal generators are the same for group 577(4), 5C/(2,2), 
and 5C/(3,1). 
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11 Discussion 

Let us summarize the main point of the present treatment. 

Parametrization of 4 x 4 matrices G of the complex linear group GL(4, C) in terms of four 
complex 4- vector parameters (k, m, n, I) is investigated. Additional restrictions separating some 
subgroups of GL(4, C) are given explicitly. In the given parametrization, the problem of in- 
verting any 4x4 matrix G is solved. Expression for determinant of any matrix G is found: 
det G = F(k,m,n,l). Unitarity conditions G + = G" 1 have been formulated in the form of 
non-linear cubic algebraic equations including complex conjugation. Several simplest solutions 
of these unitarity equations have been found: three 2-parametric subgroups Gx, G2, G3 - each 
of subgroups consists of two commuting Abelian unitary groups; 4-parametric unitary subgroup 
consisting of a product of a 3-parametric group isomorphic SU(2) and 1-parametric Abelian 
group. 

The Dirac basis of generators A&, being of Gell-Mann type, substantially differs from the 
basis Aj used in the literature on 5f/(4) group, formulas relating them are found - they permit 
to separate SU(3) subgroup in 577(4). Special way to list 15 Dirac generators of GL(4,C) 
can be used {A^} = {ai © flj © (otifij = K © L © M)}, which permit to factorize 5C/(4) 
transformations according to 5 = e iaQ! e^e e e* m , where two first factors commute with 
each other and are isomorphic to SU(2) group, the three last are 3-parametric groups, each of 
them consists of three Abelian commuting unitary subgroups. Besides, the structure of fifteen 
Dirac matrices permits to separate twenty 3-parametric subgroups in SU(4) isomorphic 
to SU(2); those subgroups might be used as bigger elementary blocks in constructing a general 
transformation SU(4). It is shown how one can specify the present approach for the unitary 
group SU (2,2) and 517(3,1). 

In principle, all different approaches used in the literature are closely related so that any 
result obtained within one technique may be easily translated to any other. There is no sense 
to persist in exploiting only one representation, thinking that it is much better than all others. 
Success should lie in combining different techniques. For instance, Euler angles-based approach 
provides us with the group elements in the separated variables form, which may be of a supreme 
importance at calculating matrix elements of the group. In turn, a factorized subgroup- based 
structure is of special interest in the particle physics and gauge theory of fundamental interaction. 
Geometrical properties of the groups, their global structure, differences between orthogonal 
groups and their double covering, and so on, seem to be most easily understood in terms of 
bilinear functions in space of linear parameters: G = XjAj, x"j = ejkix' k xi. 

We have no ground to think that only exponential functions e* A are suitable for exploration 
into group structures. We may expect that in addition to Euler angles many other curvilinear 
coordinates might be of value for studying of the group structure. For instance, in the case of 
the group 50(4, C) we have known 34 such coordinate systems owing to Olevskiy investigation 
|58| on 3-orthogonal coordinates in real Lobachevski space. 

In conclusion, several words about possible application areas of the obtained results. The 
main argument in favor of constructing the theory of unitary groups SU (4) (and related to it) 
in terms of Dirac matrices is the role of spinor methods being widely adopted in physics. Let us 
mention several problems most attractive for authors: 

SU (2, 2) and conformal symmetry, massless particles; 
classical Yang-Mills equations and gauge fields; 
geometric phases for multi- level quantum systems; 
composite structure of quarks and leptons; 
SU (4) gauge models. 
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In particular, description of the group SU {2,2) in terms of matrices a J , fP should be of 
great benefit in investigation of conformal symmetry in massless particles theory. For instance, 
classical Maxwell equations in a medium can be presented in 4-dimensional complex matrix form 
with the use of two sets of matrices, exploited above: 



dx° dxi 



M + 



* dx° ^ dxi ) 



1 

CO 



p 

j/c 



where 



M 
N -- 





M 



N 



M = —{D + iH/c) + (E + icB), 
eo 

AT = — (D - iH/c) -{E- icB). 
eo 
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